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Presentation 

Senores miembros del jurado. 

Senor decano de la facultad de ingenieria. 

Senores docentes de la escuela de ingenieria mecanica. 


De conformidad a lo contemplado por la ley universitaria 30220, Art. 45. 
Inc. 45.2, en concordancia con lo dispuesto en el Art. 133 de los Estatutos de 
la Universidad Nacional de Trujillo, presento a vuestra disposicion, bajo la 
modalidad de elaboracion de una TESIS, el presente trabajo de investigation 


titulado: 

“Esfuerzos en cab 



dos transversalmente por una 


Ingeniero Mecanico. 


como objetivo principal modelar 


Mediante el cual po 
El presente trabajo 

la respuesta de cables ;tados transversalmente por una 

masa puntual. Se ha lograc'.-' orofundiz^^’ampliar el conocimiento sobre este 
fenomeno y, en consecuencia, se ha contribuido con herramientas analiticas 
para el desarrollo y diseno de barreras de protection viales, barreras de 
prevention de caida de rocas y bases para detener objetos que caen en 
minas ya que estos disenos sirven para disminuir la gran cantidad de 
accidentes mortales que pasan cada ano en nuestro pais y en todo el planeta. 


Tengo a bien poner a disposicion este trabajo, a fin que pueda motivar la 
realization de futuras investigaciones en el tema del impacto transversal de 
cables. 
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Abstract 


A deeper understanding in the transverse impact of wire ropes is 
required to be able to design more and improved structures to reduce the 
soaring loss of human lifes in traffic accident, by fall of object in mines and 
the huge amount of rockfalls and landslides ubiquitous all along the rugged 
geography that we have. By all these reasons, in this thesis are given the 
mathematical models of the forces in the ends of a steel wire rope and the 
movement of the body that impact transversely in its middle point. To achieve 
these objectives, both ends of the wire rope are considered rigidly fixed, 
reducing its geometry to an unidimensional elastic string and modeling the 
body as a material point that a rock, automobile or any machinery 

that could fall over a nieving this way a high level of 

acceleration what in Jffi 'lacements and internal stresses 


around the contact z 
using the developed 
ropes type 3x7 and 



Next the models are validated 
5 the reaction of the steel wire 
conditions imposed in previous 


investigations done abou L orotection barji&rs. At the end is concluded that all 
models, excepting the modef^il 1 ' internal stresses, express the physics 
present in the referred experiments and simulations thus contributing to 
augment our knowledge about the transversal impact of steel wire ropes. And 
additionally, it has been contributed with a method to determine the frontier 
conditions in a transversely impacted string that can also be used as basis to 
solve more general problems, as the oblique impact of wire ropes. 


Keywords: Transversal impact, steel wire rope, mathematical model, 
protection barrier. 
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Figure 45: Initial position of the parallelepiped and base made of wire ropes an 
instant before the impact simulated by P. Gospodarczyk. The base 
has diameter of 07.5m and is used in deep shaft mines. 1 is the 
parallelepided, 2 is the set of wire ropes and 3 is a sheet metal [33], 
.89 


Figure 46: First idealization of the simulation of an impacted artificial base 
used in deep shaft mines and simulated by Gospodarczyk.90 



Figure 47: Second idealization of the simulation of an impacted artificial base 
used in deep shaftjn^i^^^d simulated by Gospodarczyk. L, is the 
distance bel t coincides with the diameter of 

the base, M Jj m a& s ^ of -lock, V00, is the impact velocity 
of the bloct d axis and O, is the origin of 

coordinates mJdjK the wire rope.90 


Figure 48: Ca 

artificial 


impacted wire ropes in an 
shaft mines, simulated by 


Gospodarczyk. Tension indicated with an arrow is the breaking 
strength.92 


Figure 49: Calculated velocity with the final velocity, from the simulation, 
reached by the 18 TN parallelepiped after passing through tha 
artificial base simulated by Gospodarczyk. Velocity indicated with 
and arrow is the calculated final velocity taking in count an standard 
breaking strength of 279 KN for a 1x7 wire rope.93 

Figure 50: Final result of the simulation made by Gospodarczyk with the block 
decelerated from 75 m/s to 48 m/s [33].95 


Figure 51: Cylindrical impact front used in the experimental car of C. Stolle 
[31].97 


Figure 52: Results of the model of the internal stresses (maximums) around 
the impacted zone of the wire rope used in the experiment 1 of C. 
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Stolle. With two curves one using the maximum calculated tension 
Tmod. Mat.Expl. Max , and the other with the measured tension 
TExpl. Max and compared with the wire rope breaking strength 
oR. Max.99 

Figure 53: Results of the model of the internal stresses (maximums) around 
the impacted zone of the wire rope used in the experiment 2 of C. 
Stolle. With two curves one using the maximum calculated tension 
Tmod. Mat.Expl. Max , and the other with the measured tension 
TExpl. Max and compared with the wire rope breaking strength 
oR. Max.101 


Figure 54: Permanent deformation of the wire rope around the impact zone 
after the Exp. 2 oJvC. Stolle i3'!1.102 

Figure 55: Permanent Vp rope around the impact zone 

.102 


after the Exp 


Figure 56: Results of ^ 



lal stresses (maximums) around 
3 used in the experiment 4 of C. 
frig the maximum calculated tension 


the impactec 


Stolle. With two%<ir/es one,u$ 

Tmod. Mat.Expl. Max , and the other with the measured tension 
TExpl. Max and compared with the wire rope breaking strength 
oR. Max.105 


Figure 57: final condition of the 3x7 wire rope after being impacted in the 
experiment 4 of C. Stolle. [31].107 

Figure 58: Results of the model of the internal stresses (maximums) in the 
impacted zone of the wire ropes of the artificial base for deep shaft 
mines simulated by Gospodarczyk using two maximum calculated 
tensions Tmod. Mat.Expl. Max , one per each side of the 
parallelepiped that impact and compared with the wire rope 
breaking strength oR. Max.109 

Figure 59: Approximated model of an impacted wire rope with and idealized 
curved line and its approximation with a triangular shape during its 
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movement and before any waves reflects at its rigid ends. M, body 
mass; V00, its initial velocity; V, it’s actual velocity; O, origin of axis 
x, y; L, is the distance between anchorage; F, represents the 
reaction force; CM, is the longitudinal stress wave speed, X*, is the 
transverse wave speed and a is the angle between axis x and the 
triangular profile sides.125 

Figure 60: Comparison of the mathematical model of the anchorages tension 
T(KN) applied to the transversely impacted wire rope used in the 
experiment 1 of C. Stolle, jointly with the simulation of C. Stolle and 
other numerical codes that he used to compare his results. [31]. 130 
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List of symbols 


A 

Transverse area of a bar or combined area of all wires in a wire 

rope, metallic area 

a 

Radio of the contact area between to cylinders with 

perpendicular axis. 

b 

Width of the rectangular contact zone formed between a plane 

rigid solid and a cylinder 

B, Bj B 2 ... 

Dimensionless constants 

c 

Sound velocity or velocity of the longitudinal stress wave in a bar 

Cl 

Transversean unidimensional deformable 

c„ 

Longitudi , elocity in an unidimensional 

Di y D 2 

Diameter^ ■ ■ > cylii^dOi ^ ■■ ■ i^^vst 

E 

Elastic modulus 

e d 

Unloading stage elastic modulus of a transversely impacted wire 

rope 


E s {Slower> Supper) Steel wire r0 P e elastic modulus equal to the slope of the 
secant between an initial and final stress 


E s ( 0, CTupper) 

Steel wire rope elastic modulus equal to the slope of a secant 

between an initial zero and a final stress 

E t,up(°) 

Steel wire rope elastic modulus equal to the tangent slope in 

some point of the stress curve while the wire rope is being 

loaded (increasing tension) 
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t,down ( 

Steel wire rope elastic modulus equal to the tangent slope in 

some point of the stress curve while the wire rope is being 

unloaded (increasing tension) 

f() 

Function 

F 

Reaction or contact force between the impact body and the 

impacted wire rope 

Fc 

Loading force along a wire rope transversely impacted 

F D 

Unloading force along a wire rope transversely impacted 

a{ 0 

Function of time 

h\ y h 2 

Curvature tn S^culate the equivalent radio of two 

cylinders^ 

i 

Natural n (j .:er|u|A 

J( 0 

Function ^ 

kv 

Transverse 'ave velocity irj,*^n impacted ideal string 

L 

Wire rope length or characteristic distance on a deformable solid 

m 

A constant coefficient 

mm' , nn' 

Cross sections of a bar 

M c 

Impact body mass 

MCable • 

Wire rope total mass 

M 

Combined mass of the impact body and of the wire rope 

n 

Natural number different of zero 

IN 

Set of the natural numbers 

n 

Contact length along a wire rope impacted by a solid with a 

plane front 
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o 


P 0 


R 


e 


^ 1 , ^ 2 , 

R, R 3 R 4 R 78 


s,S 

At 

At 

T a(0 j T b(0 


T 

± C. Max 


Origin of a coordinated system. 

Maximum contact pressure between two deformable solids in 
contact 

Equivalent radio for two cylinders in contact with perpendicular 
axis 

Rigid ends of an ideal string 
etc Dimensionless constants 



Curvilinear coordinate of reference or Lagrangian 


Actual Curvilinear coordinate of an uni-dimensional string. 
Time i 


Finite 

Times 


Maximum 


during the loading stage at the ideal 


string ends 

T D.Mdx Maximum tension reached in an loading process symmetric to 

the unloading process in a transversely impacted wire rope 

T c Loading tension at the ends of a transversely impacted string or 

wire rope 

T d Unloading tension at the ends of a transversely impacted string 

or wire rope 

T Tension in points of a string 

||T|| Module of the tension vector in a string 

tc.MAx Time taken to load a transversely impacted wire rope up to a 

maximum tension 
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t d.mAx 

Time taken to unload a transversely impacted wire rope from a 

reached maximum tension 

t„ 

Loading or unloading time of an ideal longitudinal impacted 

string without internal friction and perfectly elastic 

Ir 

Time to reflection or the time taken by a longitudinal wave front 

to come back to the impact point in a string 

^ total 

Total time taken to load (decelerate) and unload (acceleration 

and return) a wire rope after being transversely impacted 

t c 

Characteristic time of a collision, up to the reaction force 

t 

becomes maximumjDr it backs to a minimum 

u 

Displacer 

u h u 2 

Movemer IttHlj^iO^ eters of the delimiter string with 

V 00 

small disph. sis 

Impact body 'vfe/€]|fjity 5 ^^velocity at the instant of contact with 
the string or wire rope 

V 

Impacted string’s point velocity 

w(q) 

Function of a variable q 

x , y , z 

Coordinates of reference 

X 

Actual coordinate or force per unit length along the x axis in an 

ideal string 

( X,Y,Z) 

External forces per unit of length over an ideal string 

y 0 

Displacement along axis y of a string’s middle or impact point 

y R. Corr. D 

Displacement while unloading of the middle point of a string with 

small displacements 
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y r. corr. s(ed) Impact point displacement during the unloading expressed as a 
symmetric function of the loading of the same string with its 
unloading modulus of elasticity E D and with the model of small 
displacements 


y R Corr S(£D Mdx) Impact point total displacement for the string loading stage. 

yR.aprox. Impact point approximate displacement using the model of small 
displacements 

y r. corr. Corrected displacement of the string’s impacted point with a 
delimiter movement and small displacement 


y M 


y m.c 


y m.d 


y" M 

y" M . c 


Material point displacement (model of large displacements ) 


Material 
displacer 


Material 
displacer 



n hile loading (model of large 


hile unloading (model of large 


Material poim&^eleration (model of large displacements ) 


Impact body acceleration while the wire rope is loaded (model of 
large displacements) 


y '' M D Impact body acceleration while the wire rope is unloaded (model 
of large displacements) 


y ' M Impact body velocity (model of large displacements) 


y ' M c Impact body velocity while the wire rope is loaded (model of 

large displacements) 

y ' M D Impact body velocity while the wire rope is unloaded (model of 

large displacements) 

y R Delimiter string’s middle point displacement (model of small 

displacements) 
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( dyldx) 2 

Bar’s neutral axis slope value squared 

a 

Angle between the x axis and the impacted string taking a 

triangular shape after contacted by a material point 

8 

Finite displacement 

8A 

Infinitesimal change of area 

8X 

Infinitesimal change of the actual coordinate of a deformable 

solid or ideal string 

£ 

Small lineal deformation 

£ 

Deformation ratio 

tPcorr. 

Correction of^e^'^roxi mated displacement (model of 

little disp|$ %J er blue light’s middle point 

Jt 

Constant 

Po 

lnit ' al Nne 

P 

Actual lineaiVctensity of a str#ff 

Pmat 

Density [Kg/m A 3] 

o 

Stress 

O l 

Longitudinal stress beard equally by every wire in a wire rope 

O'** 

Stresses parallel to the x axis along the z axis 


Stresses parallel to the y axis along the z axis 

0* 

Stresses parallel to the z axis along the z axis. 

^l.mAx 

Maximum longitudinal breaking strength of a wire rope 

° SR .MAX 

Maximum breaking strength of a wire rope 
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o v M Equivalent Von Mises stress along the z axis with origin at the 

center of the contact zone between the impact body and the 
wire rope 

°R,Mdx Breaking strength (Pa) 

Shear modulus of two steel cylinders in contact 
v Poisson ratio 
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1 . Introduction. 

1.1 . Background and motivation. 

Protection barriers are used all around the world to reduce the fatal 
consequences of accidents in roads, by fall of objects in mines and 
landslides [1], [2], [3], [4], The reasons to design and install these structures 
are showed in the global annual statistics that show how more than 1.25 mill. 
Persons die in traffic collisions, that mining is responsible for 8% of all labour 
accidents and landslides caused 4626 deaths only between 2015 and 2016 
[5], [6], [7], [8], [9], Peru, is not an exception, to the contrary, we have much 
more reasons to use protection barriers as 50% is very susceptible to 
landslides and rockfalls (^jfA£ll5/^0cities up to 38.9 m/s) caused by 



logical phenomena (El Nino), that 
same biennium. We should also 


earthquakes and the I ;i ! 


caused 591 landslide^ A 


take into account the 


i accidents: 2253 deaths in 2016 
cal statistics in mining with 115 


(with velocities up to 


accidents from falling H]> [12], [13], [14] [15], [16]. For 



during the summer of 2017 when El Nino caused landslides and the barriers 
stopped medium size rocks. Also, in Lima, along the cliffs of the Costa 
Verde, is installed a geo-mesh to protect from rock falls to the drivers in that 
road [17], [18]. 

In figures 1, 2, 3 y 4 are showed examples of protection barriers in 
some of its applications [19], [20], [21], [22], 
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Figure 1: Example of a Figure 2: Road barrier class CASS™redirecting 
vial barrier made of steel an impacting vehicle [20]. 
wire ropes. [19]. 




Figure 3: Dynamic barrier used at ravine 
Carosio in Chosica, Lima. This stopped 
detritus during the landslides of 2017 
summer [21], 


Figure 4: Example of a barrier stopping a falling rock 
that could impact cars if it reach a road down hill 
[ 22 ], 


2 
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However all these actions are not enough to solve our situation 
because We should not limit ourselves to wait and use these solutions, We 
must do better and new designs. But to design and improve, the behavior of 
the barrier’s structural elements as the wire ropes have to be studied as 
these elements are the principal part of many types of safety barriers as: 
flexible road barriers, rockfall barriers, and artificial bottoms used at deep 
shaft mines. In all these designs the wire ropes are impacted obliquely and 
transversely by vehicles, rocks o machinery which in turn are redirected or 
stopped to avoid or lessen economic and mortal damages [1], [2], [3], [4], 
[23],[24], [25], [26], [27], [28],[29], [30], Due to this reality, this job, is centered 
in studvina the transverse impact.in steel wire ropes so, this way, We can be 



ntioned above. 


is based on investigations like 
t 1 simulated with precision a 3x7 


wire rope under transv^se impacTto bemused in the design of flexible road 
barriers, but his numerical stBijSjl^i.miied to this wire rope and low velocities 
(10.57 m/s) [31], About rockfall protection barriers, Koo et al, studied 
numerically two types of this structures to predict its behavior when a 4000 
Kg sphere impacts it at 36.6 m/s, they validated its study with experimental 
data although some simulated data differed from the measured on wire ropes 
and columns [32], Similar works were done related to artificial bottoms for 
deep shaft mines, for example, Tytko et al and P. Gospodarczyk simulated a 
design with 1x7 wire ropes in a criss-cross pattern to form a net which was 
supposed to stop machinery falling from 800m high. They were limited by the 
little knowledge about this topic and the resources to make experimentation 
at scale 1:1, hence, first they decided to test the elements or materials under 
impact individually then all this data was used to simulate and optimize their 
design [29], [33], And the more recent work about structures with wire ropes 
under transverse impact is the one of Dqbek et al, their principal objective 
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was also to design a barrier for deep shaft mines but with a fall of 1000m and 
using 6x36 wire ropes in a criss-cross pattern testing in advance its 
properties individually in laboratory. Its main conclusion was that wire ropes 
only absorb a little portion of the energy and do not optimize its design [4], 


At this point, what can be observed in the referenced researches is how 
are limited to the designs and simulated wire ropes using only a commercial 
software LS-DYNA. Neither using any analytical model nor doing a reference 
to theoretical developments about wire ropes under impact that were being 
made by many years as G. W. Housner did in 1945 when he had to model 
the transverse impact of airplanes over wire ropes attached to balloons which 
were used as a barrier to prgi^ 1 !!^^ installation during WWII, His model 

fa 

a constant velocity) could break 


predicted that the airp | 
the wire ropes with# 
practice [34], During 
theoretically the sarr 
transverse and oblige 



;e and was what happened in 
1 A. Rakhmatulin, also studied 
rage balloons, and hence, the 
is including friction, keeping the 


impact body’s velocity constant or with a small deceleration and assuming 

"6|£>r,s 

the wire rope as a string with no bending resistance and small displacements 

[35] , Then, other researches made similar formulations: F. O Ringleb, 
modeled and tested experimentally 6x19 wire ropes with a hemp core and 
used this to design aircraft arresting gear systems used in aircraft carriers 

[36] , But was N. Cristescu, who give more contributions about the impact in 
general applied to fibers and considers a broad bibliography about this field; 
furthermore he presents a numerical method to solve the case of the 
transverse impact over a string with a fixed end but assuming beforehand a 
frontier condition with a sinusoidal function which models the movement of 
the impact body (material point) [37], Then, more recently in the nineties, J. 
Lombard, based its work on G. W. Housner’s and investigated accidents 
involving military aircraft accidents impacting transversely wire ropes of 
ropeways. He establish a purely plastic deformation and assumes a parabolic 
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profile for the string instead of the triangular one used before and also 
assumes the body velocity as constant; his results coincided with what G. W. 
Housner found and either did any contrast with experimental data [38], [39], 
Then, what characterizes these works are the simplifications made in order to 
found an analytical solution after reducing the general differential equations 
by assuming a constant velocity, or giving a function for the kinematics of the 
for the impact body, or establishing a shape for the impacted wire rope as is 
shown in figure 5. 


> 



Figure 5: I 

during the i 
the velocity 
point and c , is 


a triangular shape 
velocity v , kv , is 
from the impact 
inal wave [34], 

More than 60 years hadjf§ E ji<y by to have a formulation where is not 


assumed a shape for the string or a constant velocity for the impact body. 


That was X. Teng (2005), who modeled fracture on strings like bodies 


transversely impacted by a parallelepiped, assuming perfectly plastic 


materials and considering a high impact velocity (~>500 m/s)with which 


simplifies the general dynamic equations of the string to the conservation of 
energy and conservation of momentum and founded a solution but, again, do 


not gives any experimental validation [40], 


This way in our present We have that every experimental o numerical 
study about protection barriers containing steel wire ropes limits its findings 
to the types of wire ropes used and do not use any theoretical approach 
made before, and these ones are limited to keep the impact body’s velocity 
constant, assume an small deceleration,suppose a function for its movement, 
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or even assume a priori a shape for the impacted wire rope and constant 
properties for its constitutive material (steel). However, all this simplifications 
do not happen in real applications because the principal function of a 
protection barrier is to stop some vehicle in an accident, or a falling rock o 
machinery, and thus, preventing them from reaching people in roads, towns, 
workplaces, etc. But, this implies a great deceleration of the impact body and 
the barrier protection structure should keep its structural integrity after a 
collision(s) according to its design. 



Hence, there exists the need to design more and better protection 
barriers. And knowing the limitations that simulations, experiments and 
theoretical approaches done about impacted wire ropes, this work, is 


centered in modeling 
any mass that could 
way it will be possible 
result of the accelera\ 
which do not break or 
protection barrier which is protect 


nsverse impact of wire ropes by 
ne machinery or a rock and this 
mging longitudinal stresses as a 
y assuming an elastic wire rope 
;e with the principal function of a 


Now, accordingly to the aforementioned, a research problem is stated: 


1.3 . Problem statement. 

Develop a mathematical model to determine the emerging stresses in a 
steel wire rope transversely impacted by a punctual mass. 

1.4 . Justifications. 

Justifications to work on the stated problem are: 


1.4.1 . Scientific justification. 

Because it is necessary to understand better and more the 
phenomenon of transverse impact of wire ropes. As an open field to research 
and as analytic alternative for nowadays numerical models. 
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1.4.2 . Technological justification. 

The need to predict the stresses in a steel wire rope transversely 
impacted by a mass that could be a vehicle, rock or machinery hitting some 
road barrier, rockfall barrier or artificial bottom mine, and because an 
analytical tool is needed for prediction and selection of wire ropes to be used 
in this barriers. 

1.4.3 . Social and economic justifications. 

There is also a need for more and better alternatives to diminish or 



due to landslides, rockfalls, road 
ips. Here is also included the 
lal and human resources can be 


I be done with prototypes of wire 


The prime goal is to develop a mathematical model for the longitudinal 
stresses, another of the the kinematics of the body impacting in its middle 
point and a last one of the internal stresses around the contact zone of a 
transversely impacted wire rope. 

1.5.2 . Specific objectives. 

After establishing and solving the mathematical models, it is sought to 
do the following: 

• Analyze the impact body’s acceleration dependence on the wire rope 
characteristics and its impact velocity. 
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• Study the relation of the wire rope’s points displacements behavior 
during the impact with its properties and the impact body’s mass and 
velocity. 

• Analyze the wire rope longitudinal stresses during the loading process 
(impact body deceleration) or unloading process (impact body return 
acceleration) 

• Validate the longitudinal stresses model by comparing its results with 
the experimental data of C. Stolle’s 1 research about 3x7 wire ropes 
impacted transversely and also with the numerical results of Tytko et al 
and P. Gospodarczyk’s research [291, [31], [33], 




mass. This research is limited to the impact in the point of the wire rope, 
which has an elastic behavior, with its geometry idealized as an 
unidimensional deformable string without bending strength but including a 
permanent deformation after loading to consider hysteresis. The impact body 
is represented as a material point, this is,no other effect is considered more 
than its acceleration on a straight trajectory in contact with a non-friction wire 
rope surface. This analytic approach comprises formulation of models 
oriented to predict the measures taken in experiments and simulations, i.e. a 
model of forces in the wire rope ends, a model of the impact body kinematics 
and one of the internal stresses around the contact zone. Then, validation is 
done for each model idealizing the experiments a simulations of 3 researches 

1 C. Stolle usa la nomenclatura DBC-1,2,3... para denominar a sus experimentos con cables 
impactados transversalmente y DTC-1,2,3.. para los experimentos de impacto longitudinal [31]. 
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where 1x7 and 3x7 wire rope types were used under transverse impact 
showing a good agreement, except for the maximum internal stresses around 
the contact zone. To end this work conclusions, recommendations and 
annexes are given. 


1.7 . Outline of the thesis. 

This thesis consist in eight chapters which arrangement is as follows: 


Chapter 1 critically and succinctly reviews the social problematic 
involved and previous studies done about protection barriers. There is also 
the problem statement, justifications, objectives, the scope and a thesis 
outline. 


Chapter 2 providll < 

impact of solids, dis A 

Mm, 

stress waves, dynami 
axial stresses and a 
strings. 



including basic concepts about 
-mations in strings, longitudinal 
ihavior of steel wire ropes under 
of transverse impact of elastic 


Then chapter 3, establishing- materials and methods to be used and 
followed, also gives which variables are considered along with geometric and 
physical restrictions and simplifications. Formulation and validation 
procedures are given as well. 


Chapter 4, developes the matehmatical models. 


In chapter 5, the models are validated using experimental and 
simulated data from three previous works about protection barriers which use 
3x7 and 1x7 steel wire ropes. Next, chapter VI concludes the thesis and 
gives recommendations. In Chapter VII are the references and finally chapter 
VIII includes two annexes. 
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2 . Theoretical background. 

Relevant topics for this work are: impact of deformable solids; 
deformation, displacements, stress waves and dynamics of bars and strings 
under external stresses, also,the behavior and properties of steel wire ropes. 


2.1 . ^What is considered as impact? 

Before studying wire ropes under transverse impact is necessary some 
approaches about the impact phenomena in general, as it is not common in 
an undergraduate level. 


Impact between deformable solids is observed when discontinuities in 
properties as density or oilier thermodynamic quantity are present due to 

exchange of momentum and that 
these discontinues changes can 


large changes in pres^p' 
then are propagated t 


be produced by simpl. 
as abrupt and violent 



perceive them with our senses 
>t an absolute definition of what 
terms of time; thus, any concept 


“discontinuity” or “sud ^ 

would be something subj^jva_ F^^d'kample, in the intellectual world of 
impact physics, a common discontinuity would be a wave front traveling from 
the contact zone to every point inside a solid that could be a bar colliding with 
another one, provoking high peaks of pressure in just some tenths of a 
nanosecond. Another example can be observed in geology, where 
earthquakes can produce severe and non-reversible changes in many 
seconds. Curiously, these cases, even happening with time and space scales 
too different, deal with the same orders of magnitude when are related with 
its life spans which are taken as its characteristics magnitudes: suppose that 
a solid medium has a characteristic space magnitud, L , (a bar length for 
example) and a characteristic elastic stress wave velocity of propagation, or 
sound velocity. If this body is subjected to an external dynamic load with a 
characteristic time t c which is, for example, the time needed to reach its 
highest level or the time when its value becomes null again; and if t c «:L/c , 
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possible to determine 
structural changes in 
criterion is considered 


then, the stress distribution and deformation in that solid are not uniform. 
Hence, wave propagation effects should be considered. Some examples of 
characteristic times are 10 nanoseconds that is a 4.2 x 10' 18 part of a human 
life of 75 years (or it could be a generation or a structure’s lifespan) while 1 
second could be a characteristic time, as it is a 6.9 x 10" 18 part of earth’s life 
(or can be compared with a geological like Cretaceous or Mezosoic). These 
very little times are also why many models about solids impacting strings 
assume that after the first contact the string points change its velocity 
instantaneously from zero to the body’s velocity and is represented with a 
Heaviside function [36], [40], [41]. However, all the above only tells that every 
initial contact between solids generates a propagating stress wave and if its 
characteristic time is ' ^t is^gaisidered a percussion 2 . But is not 

ig enough to cause considerable 
iolids. For this reason another 
■acterization”, there is suggested 
:act between solids is clasified, 


some references to c. term in; 
based on its deformation ratios ana relative impact velocity. If the observed 
deformation ratio has 10° ordehiisfis^rgnitude and a velocity of 50 m/s, then 
only a localized plastic behavior is could happen, but with velocities from 50 
to 500 m/s and deformation ratios with an order of magnitude of 10 2 then only 
a plastic behavior is expected. Of course these parameters are only 
referential as them do not consider many other factors like geometry, 
generated heat, friction or strength similarity between the colliding solids 
materials and its behavior with high deformation ratios as happen when a 
hard steel impacts an 1100 aluminum, this latter will develop a purely plastic 
flow typical of hyper-velocity impacts while the steel will have at most, an 
elastic disturbance [42], [43], Hence, it is imperative to analyze every case 
and to establish which parameters are relevant. An example is how N. 
Cristescu, pointed out that if a nylon strings under transverse impact is hit by 
a body with a velocity exceeding 1/20 of the transverse wave velocity or 



2 Synonym of impact or collision, at least in mechanics. 
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deformation ratios of 33.33 m/m/s then the effect of this wave propagation or 
the local inertia forces must be taken into account. This is to say, here one 
can not simplify the string movement equations with a priori suppositions like 
assuming a triangular shape during the movement as other have done [37], 
[34], [38], Although until now there is not a similar criterion for steel wire 
ropes that can tell when local inertia is relevant according to a impact 
velocity, this remains relatively unexplored. 


2.2 . Strings displacements and deformations. 

It is also necessary to clarify the displacements and deformation 
concepts relative to this kind of structural elements because a wire rope can 


have large displacement^ 
Love indicates how 


small displacements 
can be small enough 
Strength of Materials t 





s f p r mations. About this fact A. E. H. 

% 

on strings, can have relative not 
even with this its deformations 
its of a mathematical theory in 
Thus, there is an approximation 


using “small deformat ,obs which, lets consider a linear elastic response for a 
material and in the otr'Sfflj^g^y liiere is the “small displacements” 
approximation where elemental volume rotations can be considered as rigid 
body rotations and be irrelevant to calculate elastic deformations, but if taken 
into account they introduce non-linearities in our equations [45]. This latter 
condition can usually be found in conventional engineering applications, the 
classic example is the simplification of the bending moment - curvature 
relation after being linearized by considering the square of its slope ( dy/dxf 
too small, compared to 1, t but this can only be done after considering small 
deflections and an almost straight beam or bar [46], Thus, in any other case 
involving slender bars or wires where its deflections are comparable with its 
own length an elemental linear theory can not be used. Hence, to model 
bars, wires and strings shall be known that large displacements will be 
present. Which means non-linearities due to its geometry but with a material 
constitutive law that can be used to restrict the study to small deformations. 
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The foregoing have to be complemented by expressingmathematically 
the deformation of an uni-dimensional body which is an string as is the elastic 
solid studied here. There are various deformation definitions to choose from, 
but for a small deformations regime (in case of elastic metals) it is possible to 
probe that every one ends to be equal to Cauchy’s definition [47], It have to 
be clear too, that even if an impact involve dynamic stresses this does not 
necessary implies the material behavior would differ from its statics 
response. Von Karman, pointed out this on its study of longitudinal impact 
reaching plastic stresses in bars. He noticed how the ratio of deformation can 
influence stresses only at very high impact velocities, but as does not specify 
how high then the Weihrauch characterization can be used to establish a 
reference [43], [48], [4 


2.3 . Unidimer 

These appear 
impact. It is only need 



iteraction is not considered an 
position between its points. 


Aiming to find a mathematical expression, figure 6, represents the 
model of prismatic bars, or straight bars, there is a bar with variable 
transverse area along its axial axis aligned with the coordinate reference axis 
x, with origin O [50], [51]. Then with an isocoric process and in dynamic 
equilibrium the following can be written: 


w, . \ * „ d * u d * f r * d * u \ ( A 0 d x u ^ 

A(l + e) = A 0 ; A 0 p maf —= — (EA —)= — (£ —— ) 


( 1 ) 


d t t d x x d x x d x x l+ed x x 
In equation (1), a change of the transverse section area is considered 
related to the deformation in the actual bar configuration, but as indicated by 
Timoshenko, this does not take into account lateral inertia from particles 
moving due to contraction of these sections, hence,this relation is valid only if 
the stress wave amplitude is large compared with the transverse bar 
dimensions [52], This equation can be done simpler by considering a solid 
bar with an uniform transverse section, as follows: 
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0fu 2 

= C 


dlu 


d t t 2 d x x 2 


; c=-\l Elp w 


Being its general solution the one of D’ Alambert: 


( 2 ) 


u(x,t) = f 1 (x-ct)+f 2 (x+ct) 


( 3 ) 


X 


x 


m 


u 


n 


, d xU 

U + ——— OX 

OyA 




A+SA)(a+^^SX) 


Figure 6: Dynam deformable bar with varibale 

transverse section, m , is the coorg^te of reference, X , is the 
actual coordinate; SX is ( .^ longitudinal infinitesimal element, mm' 
and nn ' are sections, u is tne displacement, 8 A is the infinitesimal 
area’s change, o is the evenly distributed stress in the cross section 
and c is the longitudinal stress wave velocity [50], [51]. 


2.4. General equations of motion for an elastic 
string. 

N. Cristescu was the first to work completely and developes in deep the 
general equations of motion for an extensible string moving in a tri¬ 
dimensional string [37], [53], He used a curvilinear coordinate S along the 
string with an arbitrary origin, as shown in figure 7, this is an Eulerian 
coordinate and called “actual” as it measure the deformed string in the actual 
time t. The total force in any string point is denoted with T , (it is force and 
not stress because an unidimensional string has not area), is tangent along 
the string and also varies continuously with time. There is also considered a 
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segment between s 1 and s 2 tensioned and with external forces ( X,Y,Z ), 
which are proportional to the segment length. 


J 


T(s 2 ) 





\ + X = p 


z . Here, p is the actual string 


Figure 7: 

a trid 1 


string segment in 
internal and 


as follows: 


density and function of s and t . Projections of T over the three 
coordinated axis are ||r||(d s x/d s s), ||T||(S, 1 y/S,s), ||r||(d s z/d s s) , being 
||T|| its modulus and s l , s 2 arbitrary, s can also be expressed using the 
Lagrangian or material coordinate s 0 , in the same way for x , y , z there is 
a expression like: x(s(s 0 ),r), etc. Then: 


d s x 0 v oLs 


, idem for y, z. 


d So s o d s s d So s o 
Now, deformation is defined following Cauchy’s: 


d.,s-d..s € _l 


1 Pi t 2 

d s.y\ , 

Z 

2 a 

O , X 

i . 




d s s 0 

* > 



Idem for y, z 



( 5 ) 


( 6 ) 


15 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Peru. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 


























BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DE INFORMATICA Y COMUNICACION 


And due to mass conservation: pd s s=p 0 d s s 0 => p 0 =p( 1+e) here, p 0 
, is the non-streched string density. Hence, the movement equations (4) can 
be expressed as follows: 


d.x . . 

-p 0 —^+(l+e)X=0 , idem for y, z. (7) 

d t t 

Here, a material’s constitutive law can be added. This is usually an 
increasing monotone function and defined as ||T||=||T||(e) . Then using (6) 
and equations for x , y and z from (7) the next system is found: 


d s s 


llrll 8 - x 

1+e d s s ol 


d s x d] x d s y d] y d s z d] z 
ds 0 *od s / 0 d So s 0 
dlx 


(1+e) d s 0 s o 

With similar equ 
equations is hyperb 
characteristic curves, 



/ x d s e 
(l-c)^- = 0 


s 


o 


'-H^l + ( 1+e ) X = 0 


( 8 ) 


3, s 


0 


rhis partial differential system of 
d has 4 different families of 


d t s„ 

iirii r 


[1 d F \\T\\] 

d t t~~ 

.Poii+e) 2 ] 

’ d t t ~ 

Po d £ e 


1/2 


( 9 ) 


This is why characteristic lines are generally curves with a slope as a 
function of ||T|| and e . From a mechanical point of view this lines represent 
wave fronts, then there are two propagating waves with the following front 
velocities. 




p 0 {l+e) 2 


1/2 




1 d. 


Po d„e 


1/2 


( 10 ) 


Here c 7 is the transverse wave’s front velocity and only change the 
string shape but does not deforms it, that is, a material point reached by a 
transverse wave acquire only a velocity to its tangent. But if a point is 
reached by a longitudinal wave front with velocity c u it cause only stretching 
but not modify its slope. 
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2.5 . Modulus of elasticity of steel wire ropes 

As the foregoing theory is used to model steel wire ropes, it is 
necessary to know how this mechanical elements behave. Elongation in a 
wire rope depends, of course, on how its parts are arranged and what 
materials are used (Figure 8) specially its wires, because, if its core is 
synthetic then it is not taken in count as does not contribute greatly to its 
elongation properties. 



Central 

Wke_ 


Figure 8: Steel wire rope 
elements [54], 


Generally speaking, the elasticity modulus of a wire rope is very 
different from its materials modulus. They stress 3 -strain curve is usually non¬ 
linear, this is, given a type of wire rope its modulus would not be constant 
above certain stress and is a function of the level of tension. This is why in 
practice is used an average elasticity modulus. Although, simple spiral ropes 
made only by twisted wires have a minimum non-linearity and can be 
neglected in most cases [54] [55], For simple spirals and some wire ropes 


3 Stresses in a wire rope are defined as o=T/A, here A is the metallic area or a combination of 
transverse areas of every wire. 
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made of many strands with rounded wires it is possible to approximate its 
modulus of elasticity using analytic methods only, but this is not too exact for 
most wire ropes with more elaborated designs and experimental data is 
needed in order to approximate its modulus that depends on what level of 
load and how was loaded before. Additionally, the elastic limit can be defined 
as that load from where a stabilized wire rope cannot reach its previous 
length after unloading. Then, some definitions of the elastic modulus which 
have practical relevance are [54] [55]: 

• E s (a lower , cr upper ) , this is the slope of a secant line between an initial 
and final stress. It is used with fluctuating loads. An special case 



inside and also friction takes place with every loading and unloading, not only 
elastic deformation occurs. For these reasons it is recommended to stabilize 
a new wire rope, before making any measure, with 10 loading cycles. It is 
also possible to use the tangent slope in every point but in practice is only 
used as an evaluation parameter to make comparisons. Figure 9 shows the 
geometric interpretation of every modulus definition, slope of line OB is 
modulus E s (0,a upper ) between a zero stress and the maximum at which the 
wire rope works. If there is an initial stress (point E) and a final one (point F), 
then the elastic modulus would be the slope of the secant line EF and equal 
to E s (o e ,ct f ) . . Also, in the same figure, tangent lines Ot and in point A 
E t , 0 (cr),E, A (cr) . And here a correct or practical modulus would be that 


18 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Peru. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 





BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DEINFORMATICA Y COMUNICACION 


measured as the slope of line CB where the wire rope has been unloaded 
[54], [55], 





Figure 9: Elasticity modulus definitions as the slope of 
secant lines of OB, CB, and EF, and of tangents in O or A. 
cr, is the stress, e t , is the final deformation, e c , is the 
deformation while loading y e s , is the deformation while 
unloading cp , are the angles of the tangent and secant 
lines. [55]. 

Next, Figure 10 illustrates, in a general way, the non-linear behavior of 
many new steel wire ropes when being subjected to one cycle of loading and 
unloading. Plastic coated wire ropes or with a fiber core tend to have a 
completely non-linear behavior while spirals (strands) have a predominantly 
linear response. An important observation here is how after a first cycle every 
type of wire rope presents an increment of its elasticity modulus, this is why 
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10 cycles are recommended before an stable measure is taken. In the same 
figure the lineal idealization of a loading and unloading cycle is represented 
for a new wire rope using straight lines [31], [54], [56], [57], [58], [59], [60], 



Figure 10: Curves 
wire ropes, after 
[57], [58], [59], 


of some types of new steel 
first time. [31], [54], [56], 


Steel bar not drawed. 


Linear idealization of a 
steel wire rope. 


Wire ropes of types lx, 
3x. (with only 1 strand 
or 3). 


Wire ropes with more 
than 3 strands: 4x, 6x, 
8x,..18x..., filled with 
plastic or with a fiber 
core. 
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2.6. Physics of elastic strings under transverse 
impact. 

Here, only geometry and English are used, avoiding mathematics on 
purpose, to communicate the involved physics happening when a body 
transversely impacts a wire rope, and also how its quasi-isomorphic 
idealization, that is, an uni-dimensional string without bending strength 
behaves after being impacted by a material point. 

First, physics of a transversely impacted steel wire rope can be 
describe as follows: 



1. When a body impacts transversely a wire rope then the initial contact 


is produced over.fi 


this first instant 


wire ropes coi 
stresses, shear 


or/rn^^isuperficial wires which in turn will 
all along the wire rope. From 
velocity between the body’s and 
:ontact also means there surge 
Sf>c art from the longitudinal ones. 


deform and tra 


2. The longitudinal ci formation appf^ng in each wire correspond to the 


one of a double or si 


elix body (as in a spring), depending of 


what type or location it has. This means, there is a tension, torsion, 
bending and shear stresses [54], 


3. Because of the foregoing, surging stress waves are also of torsion, 
bending and also of shear nature. Every one of these perturbations 
travels through the curve axis of every wire and they add up each 
other to produce the wire rope’s total displacement. In addition, over 
the surface of every wire appear shear waves due to friction. And all 
them travels up to each wire rope end where refraction and reflection 
occurs. Then every reflected wave goes back and combines with the 
actual stresses in every wire changing continuously as the contact 
force or reaction force over the impact body changes. 
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4. Intern friction forces between wires are not a function of a constant 
coefficient only but an state function of each material in contact [61]. 

5. Also, medium resistance and gravity should be taken in count as them 
influence the wire rope geometry. 

Hence, if any research considers not to simplify the wire rope geometry, 
the model should include all deformations with its internal and external forces 
mentioned before and how reflection and refraction happens in each end. 
Interaction with the impact body should be considered too, as its dimensions 
can be comparable to the wire rope length and because friction arises in the 
contact zone. But this way is not at the author’s reach, by now, therefore here 



a wire rope is modeled a: 
material point. 


I string and the impact body as a 


/m 

Thus, in order J 

JfsM 

simplifications are mai 


rope movement these general 


1. The wire rope it \ T- : v trical axis which initial state is an 

uni-dimensional curve, straiaht with uniform density resting in a plane. 


2. Only the axial deformation is considered, choosing the Hooke’s law to 
relate it with the axial stresses. 

3. The impact body is reduced to a material point with a given velocity 
which starts to decelerate when it touches the wire rope. 

4. It is considered that the point where the body impacts, acquires the 
same velocity instantaneously after the first contact. 

This last simplifications is a consequence of assuming a deformable 
and with no bending strength string. It is a discontinuous change of velocity 
which also means the geometric model has a singular point, kink or 
discontinuity on its tangent at the impact point. This singularity is where a 
transverse wave front begins provoking that all contiguous points gradually 
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get a transverse movement which means a continuous change of the total 
mass in movement. Its velocity is represented with the symbol C t of 
equation (10). 

In the other hand, propagation of stretch is done longitudinally and its 
velocity depends on how much deformable is the string. Here depends on its 
elasticity only and represented by C n of equation (10); it separates relatively 
every point it reachs but without altering its tangent direction as the 
transverse wave does. 



Figure 11 shows how the shape of an idealized wire rope changes 
(before any reflection takes ^gj^^fter being impacted transversely by a 
punctual mass. The t • has a velocity C, which varies 

with tension, The lonf 4)° c ity depends on how tension 

changes with respec is, ' n e l as ti c case it is a 

constant C n . The irr llg t^(| p ^^Tocitv. at the instant of impact, is 


V 00 , the body has a 


actual velocity, S 0 , is the initial 
nt and S is the actual one. The 


coordinate (Lagrangian) 
space where the movement occur does not offer any resistance or friction. 


MTV 

oo Transversal wave front 


It 


Longitudinal wave front 



i Wire rope idealization as an 
V Unidimensional string 


Figure 11: Movement of a string transversely impacted by a punctual 
mass. M, is the impact body mass, Voo, is its initial velocity, V, is its actual 
velocity, O, represents the origin of the coordinated system x, y, S, is the 
actual curvilinear coordinate and So, is the referential coordinate. 
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Ideally, physics of the transverse impact of a material point with a 
deformable unidimensional string has the following stages: 

1. Stage 1: At time zero an initial contact is produced between the 
material point and an point in the string, and instantaneously both get 
the same velocity. This way the loading stage begins. 


2. Stage 2: Here the body starts to decelerate and from the contact point 
a longitudinal wave travels with the velocity of sound along the string, 
at the same time a transverse wave, that only changes the tangent of 
every point that reachs, propagates to both ends. This last wave is 


how local inertia appears although the longitudinal wave also 
contributes but 'r^pmrior - ;• 

rw '% 

% 


changes its 


3. Stage 3: If the 



^ause it stretches the string and 


lecelerates and the longitudinal 
effected waves are added to the 
ie reaction force over the impact 
body, this is, the stress wave produced initially in the contact point 
comes back and influer?c^the movement. The transverse wave also 


wave fronts re£ 


stress state the 


reflects or refracts and begins to influence the string points and body 
movement as it interacts with the wire ropes inertia. 


4. Stage 4: The wire rope gets a zero velocity and begins to change the 
sense of its movement (keeps its straight direction). This way the 
loading stage finish and the begins to unload. Then, should be 
considered if the impact body is kept attached to the wire rope and if 
this one has a permanent deformation, at some time while returning, 
the point of contact will have a minor velocity than the material point 
and its movement will continue until its energy totally dissipates. Here 
the movement is considered only until the wire rope tension becomes 
zero again indicating the end of the unloading stage. 
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3 . Methods and materials. 

3.1 . Information gathering. 

All the necessary information to establish the problem and solve the 
equations is found, mainly, in works done by F. O. Ringleb, N. Cristescu and 
K. Feyrer [36], [37], [53], [54], Validation with experimental and simulated 
data is done by using C. Stolle, Tytko et al and P. Gospodarczyk’s 
researches [31], [29], [33], 


3.2 . Object of study. 



Here is studied a wire rope anchored in both ends, made totally of steel, 
helically braided, and with circular cross section wires. It gets a movement 
after being impacted traf^erse^^^punctual mass. And in order to 

geometry and movement are 
simplifications. To validate the 
n C. Stolle, Tytko et al and P. 
ly impacted wire ropes type 3x7 


analyze mathematical^ 
idealized applying m 
established mathema 
Gospodarczyk researc 
(ASTM A741 -11 (2016 ft and iy po i a i 
[33], 


5896-2012) are used [31], [29], 


3.3. Methodology to establish the mathematical 
models. 

The method used in this research is part of the mathematical modeling 
methodologies and these in turn part of general modeling which main 
purpose is to relate the real world with some simplified physical and 
conceptual representation moving back and forth cyclically between reality 
and imagination. Specifically this work is within the Philosophy of change 
founded by Newton and developed, here, by J. Olivencia 4 . 

The first step to establish a mathematical model is to set restrictions 
and simplifications based on the real phenomenon and then mathematics 
can be used as language to express what is in the idealized world [62], In this 

4 Professor of mathematics at the Universidad Nacional de Trujillo, Peru. 
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first part, the idealization is done following to P. Teodorescu. Thus, 
observations of transverse impact in wire ropes in many types of protection 
barriers are translated to an ideal reality in the mechanics of deformable 
solids modeling. This way, any physical bodies using its more simplified 
shape: bars are reduced to its axis, plates are reduced to its middle surfaces, 
or a wire rope is taken as its longitudinal axis. Hence, all is framed into the 
Strength of Materials [63]. Then at the ideal world and needing to 
establishing a theory, the Rational Mechanics paradigm can be followed to 
set specific steps to model mathematically a wire rope impacted transversely 
by a punctual mass. So, following what A. Klarbring pointed out the method 


is: [62]: 



1 . 


alized setting geometrical and 
v cations. In this way the bodies 


hl a geometrical and mathematical 
Is for material bodies: discrete, 
£ id three-dimensional ones. 


2 . 


should be included: mass, energy, 


linear and angular momentum conservation; then, expressions in 
terms of differential equations can be written accordingly to the 
bodies under study. 

3. And the third step to obtain a solvable mathematical expression is 
to introduce particular laws: physical and mathematical 
restrictions, force and constitutive laws (that can be experimental 
or idealized) 

Furthermore, a validation methodology implies a quantitative and 
qualitative comparison between the experimental or simulated results and 
what is predicted with the developed theory. And then analyze why there are 
differences and coincidences, taking in count how the chosen parameters for 
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this model, the wire rope and impact body characteristics influence the 
analytical and experimental results. 

3.4 . Research variables. 

If the mathematical models of this project are classified as described by 
K. Velten y a F. E. Cellier, the system in study in each model is: technical, 
deterministic, continuous, bi-dimensional and physical. Knowing what has to 
be answered each model is mechanistic, instationary (in time), distributed (in 
space), inverse, of research (not management), analytical and in macro¬ 
mechanics. And the mathematics is: not-lineal, numerical, continuous, and 
with differential equations [64], [65], In sum, all models are gray-box models 



c ■ 


nt of impact, V 00 . 


3. Wire rope’ JfJOk joke’s law for the loading and 


unloading st&/p rr— f. f 



4. Wire rope Poisson ratio, v. 

5. Wire rope transverse or metallic area, A . 

6. Wire rope equivalent diameter, D . 

7. Wire rope lineal density, p . 

8. Wire rope straight length: L . 

The controllable variables are: 

1. Wire rope elastic limit. 

2. Wire rope breaking strength, a R Mdx . 
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Uncontrollable variables are:. 

1. Margin of error of the solution method. 
Output variables are: 

1. Wire rope stress or force, o, T . 


2. Acceleration of the material point, y '' 


M ■ 


3. Velocity of the material point,, y ' M . 


4. Displacement of the material point, y 


M ■ 


5. Equivalent stresses around the impact zone the wire rope, cr, 

. ;>-‘ J ' r % 
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Figure 12: Grey box representation of the established 
mathematical models. 


28 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Peru. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/Hcences/by-nc-sa/2.5/pe/ 







BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DE INFORMATICA Y COMUNICACION 


3.5 . Procedure. 

The procedure follows the outlined method before. This way, first 
geometric and physical restrictions and simplifications are established. Then 
a procedure is formulated to establish the mathematical models, its solution 
and the validation with experimental and simulated data. Physical and 
geometric simplifications and restrictions. 


This study has as reference the first impact scale from the Weihrauch 
characterization, this means plastic deformations are localized (contact 
zone), impact velocities are under 50 m/s and deformation ratios with an 
order of magnitude of 10° mm/mm/seg are expected [43]. Restrictions and 
simplifications for the stee! r ^jSr#'ltl^fetr^the following: 


1. The wire rope 
zero initial tens 
as a constant p 

2. Only wire rope 



o an unidimensional string with 
:y. A transverse area is assigned 

modeled, with helically braided 


wires and a circular section and without natural or plastic core. 

3. Only linear effects from the stress-strain curve are taken in count and 
up to the wire rope breaking strength 


4. Effects of torsion and shear stresses in the wires are not considered. 
While its friction and wires rearrangement, during the loading and 
unloading stages, are modeled using two different elasticity modulus 
for each stage. 


5. Small deformations are supposed and only dependent of time but not 
of position, temporal ratio or temperature. Its material is considered 
isotropic. 


6. It is also supposed that the applied tension does not change the wire 
rope volume, this is, the steel suffers an isochoric process and does 
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not present any shock wave that could provoke a phase change of the 
materials considered. 

7. The impact is produced in the wire rope middle point or coordinate 
zero in the string. 

For the impact body the following is used: 

1. All bodies that could impact as: rocks, vehicles or machinery are 
assumed as material points and with a known mass and velocity, 
although in the contact stresses model it is assigned a shape to allows 
calculations. 

O l+o mo*\/oiman+ o+roioih+ onH fno\/nrpn +/", +l-m iA/irQ 



velocity instantaneotfejr/fjsrifisS^’to the body’s velocity. 


An idealized wire rope impacted by a material point is presented in 
figure 13, mass of the material point is M , its initial or impact velocity is 
V 00 , the reference coordinates are S 0 , actual coordinates are S and the 
body’s actual velocity is V .Procedure to formulate the mathematical models. 


Rigid anchorages ^ : ^oo ^ 




y? 

L 


Figure 13: Geometry of an idealized wire rope transversely impacted by a 
punctual mass. M, is the impact body mass, V 00 , is its initial velocity, V, is its 
actual velocity, O, is the origin of coordinates x, y, S and S 0 , are the string’s 
actual and referential curvilinear coordinates L, is the wire rope length. 
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1. First, a continuous time t is defined, symmetric in every moment 
(Thus, conservation of energy and mass can happen) and independent 
of the space or existing bodies in it. An space is also defined from an 
inertial reference with cartesian coordinates Oxy , euclidean 
properties, bi-dimensional, continuous and symmetric (Thus, 
conservation of momentum can happen). 


2. In the defined space is a steel wire rope idealized as an uni¬ 
dimensional finite solid, deformable, with uniform density, without 
bending resistance, of elastic material, initially straight and coincident 
with x axis. 


3. 


4. 



There is also a 
idealization of tbs* 

/IF 

velocity V 00 at i 
and it is presupj 
the loading stag 


^ral^pofM^fifined with 


mass M , and is the 

ihes the wire rope with an initial 

o 

gative direction along the y axis 
arge deceleration at least during 


After defining the geometry, the general laws of mechanics are used to 
establish mathematically the general model, these are given in 2 the 
theoretical background, hence, here only the constitutive law is 
needed. This law is defined with an equation relating the stress, strain 
and the elastic modulus a=Ee or its equivalent with the string 
tension using Hooke’s law T=AEs with a modulus of elasticity if it is 
loaded and a different when it is unloaded. 


5. Next, the equations of motion are stated for the impact body, doing 
mathematical simplifications in concordance with the physical 
conditions. 


6. After finding numerical answer using the experiments and simulations 
from the references mentioned, a model is established to predict the 
stresses inside the wire rope around the contact zone, or at least the 
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maximum ones. To do this Hertz theory is used and a geometry and 
material is assigned to the impact body. 


3.5.1 . Procedure to validate the mathematical 
models. 

Here only a quantitative comparison can be made for the longitudinal 
stresses at the anchorages and the body kinematics because the 
experiments and simulations referenced only present these data. The contact 
stress model can only be compared qualitatively with what was observed 
because there are not numbers recorded for each experiment or simulation 
about the internal stresses in the wire rope around the contact zone. The 


following is considered: dP 


1. If the wire rope 
curve have an 
then there will b 



ated is new, then its stress-strain 
;tion and wires rearrangement), 
1 by a small non-linear transition 
to the plastic zo ^¥$ 1 $ before here, as was established 

before, is considSyfed only the linear part but keeping the same 
breaking strength whicn was measured (in the referenced works) or is 


standard. 


2. Calculations are done from the first contact or after the simulation or 
sensors detect a tension bigger than zero and up to it is zero again or 
the wire rope breaks. 

This is how, the longitudinal stresses model and the contact stress 
model are validated with the experimental results of a 3x7 wire rope under 
transverse impact done by C. Stolle and with the simulations of a 1x7 wire 
rope under transverse impact done by Tytko et al and P. Gospodarczyk [31], 
[29], [33], 
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4 . Mathematical modeling. 

4.1 . Modeling the longitudinal stresses. 

To solve the differential equations (6) and (7) which represents the 
general movement of an uni-dimensional string, many methods have been 
developed along the past century 5 . The last one is that of Yu. A. Dem’yanov 
(1993) who presented an asymptotic method (without giving examples or an 
experimental validation) based on the method of characteristics to solve the 
problem of wave propagation over a transversely impacted string by a 
material point [66], In order to find an answer he makes some simplifications 
like assuming the string deformation and the material point deformation as 


functions of time as follows. In figure. 13 are the variables, Bi are constants: 


t(' 

So, if there is a 
does not limit the po 
used with 



( 11 ) 


Ji>0 (12) 

with small deformations (Which 
acements) and equation (11) is 


Hooke’s law T~EAe , then, ttei^ovement equations are: 


d,s 


Ill’ll 


1+e d.s 0 


-Po 


d]x 

dj 2 


:= 0 


d<s 


Ill’ll 


1+e s n 


-Po 


$y_ 

d t t 2 


=o 


(13) 



And both equations can be solved with a change of variable [67]: 

x, y = w{q); q = \a{2+m) 2 (S 0 +B 2 Y- 4t 2+m \ m ;a = -^ (15) 

A common solution is given for both: 

5 Before, many others did experimental and analytical research about strings (wires) under transverse 
impact using other methods to solve the partial differential equations (6) and (7), but assuming a 
constant velocity for the impact body or as an exponential function of time. See: “Experimental 
evaluation of material Behavior in a wire under transverse impact”, A. B. Schultz, P. A. Tuschak, A. 
A. Vicario, ASME 1967, “Determination of the stress-strain relations in a wire under transverse 
impact” A. B. Schultz, T. C. T. Ting , ASME 1968, “Determination of the unloading boundary in 
transverse impact of an elastic-plastic string”, P. A. Tuschak, A. B. Schultz ASME 1972. 
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R = 


w=B 3 e 
-m/2 


~Br 


2+m 


B,=— 1/am 


( 16 ) 


from a solid to a matr r : 
without bending resist 


In the other ham 
validate the models, c i 


Yet this solution for the displacement of the impact body at 
> y(s 0 = 0;r)= w(s 0 =0; t) is different of what Yu. A. Dem’yanov presupposed. 
Also, if it is considered that the body has a straight trajectory, this is, 
0 So x/0 ( t(s o = O;r)=O , then, the solution predicts it has a velocity along axis 
x . Thus, here is a contradiction of what its assumed by the proposed 
method and what results by using its assumptions. And even more it does not 
give a solution about how the unloading stage has to happen or how it could 
be modeled. This incoherence has its roots in how the geometry is idealized: 

'$( 5 ® -ope to an uni-dimensional string 
^ i presented in Annex 1 6 . 

simulations, referenced here to 
the ends of the wire rope with the 
acceleration, velocity the i m P a ct body. Hence, these 

frontier conditions have %be determined-first, and with that the stress state 
is determined in every point and in its critical points: ends and impact point. 
In sum, the objective to know in deep the movement of a transversely 
impacted wire rope have to be delimited to what was measured and 
calculated experimentally, and in general to how the engineering application 
is made. Then, with these formulations proceed to make a validation. This 
means the specific searched values are the forces at the ends and the 
displacement, velocity and acceleration (or force over it) of the impact body. 
After validating these results the stress values inside the wire rope can be 
approximated considering what predicts the internal stresses model around 
the contact zone. But, first, to reach these answers it is necessary to observe 
specific and common physical characteristics of each experimental and 
simulated case [29], [31], [33]: 



6 N. Cristescu avoided this problem by assuming a priori a sinusoidal function to model the impact 
body displacement as a frontier condition and then found the string movement and stresses [37], 
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1. The mass of the impact body is very large compared with the wire 
rope mass * 1 2 * * * * 7 . 


Research of: 

Wire rope 
type 

Wire rope mass 
(Kg) : M . „ 

Impact 
body mass 
(Kg): M c 

Ratio 

M /M c 

c wireRope 

C. Stolle 

3x7 

1.06E+01 

7.76E+02 

7.34E+01 

Tytko et al. 

1x7 

1.17E+00 

4.00E+03 

3.42E+03 

P. Gospodarczyk 

1x7 

8.78E+00 

1.80E+04 

2.05E+03 


Table 1: Comparison between the wire rope mass and the 
impact body mass. 


2. Time to reflection, or time taken by the longitudinal stress wave front 
to go from the impact point to an end (anchorage) t R =0.5L/C n , is too 



alculated 
time to 
flection t F 
(Seg) 


5.00E+03 


small compared with tte^jjjigfjon of the movement lifespan or loading 

fer^l 


and unloading 
ttotai/t R in table 2^' 


Reasearch of: 


Wire 

rope 

type 


C. Stolle (exp. 
1,2-4) 


3x7 


Tytko et al. 


1x7 


relation is shown with the ratio 


Movement 
lifespan t lolal 
(Seg). 


Ratio KoJ 


1.36E-03 


3.00E-01 


2.21 E+02 


1.00E-04 


3.80E-02 3.80E+02 


P. Gospodarczyk 1x7 5.00E+03 3.75E+00 7.50E-04 1.00E-01 1.33E+02 


Table 2: Time to reflection taken by the longitudinal stress waves in a transversely 
impacted wire rope compared with the wire rope movement lifespan. 


Then, in accordance with the ratios of tables 1 y 2, it is possible to say 


that in the described movements the following happens: 


1. The wire rope inertia is too small compared with the impact body 
inertia. This means that reflection of the transverse waves does not 
has a major effect neither in the forces inside the wire rope nor the 
impact body movement. 


2. As the time to reflection is too small compared with the total time of 

loading and unloading, Wai, of the wire rope. And apparently at first 


7 For P. Gospodarczyk’s simulation the mass of every wire rope in its simulated base is 

considered, this is 420 1x7 wire ropes of 7.5m length and 0l5.7mm each one. 
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glance it could be considered that every time a time, t R , passes, the 
stress waves make the wire rope tension uniform along it from each 
end to the impact point. 

Hence, point 1 permits to establish the equations considering a 
concentrated mass of the wire rope at the impact point instead of using a 
distributed one, because its inertia does not contribute significantly to its 
movement. Point 2 says that the time taken to have almost uniform tensions 
along the string is too small that the tension in each end is very similar to the 
tension in the impact point. However, because the wire rope shape, during 
the movement, is not a perfect triangle, then, the force at the impact point 


has a different direction depending of what angle has the wire rope there. 



Figure 14: Model of delimiter ideal strings delimiting the idealization of a 
transversely impacted wire rope on its central point. M, body mass, Voo, its 
initial velocity, V, it’s actual velocity, O, origin of axis x, y, S 0 , reference 
lagrangian coordinate R ly R 2 are the anchorages. Y R , displacement of the light 
blue line. Y M , displacement of the green string’s middle point (or the material 
point) and L is the distance between anchorages. 

In figure 14 is used an analog concept to that of frontier conditions but 
here are presented geometrically as “delimiter movements”. This way the 
string movement is delimited between the movement of two strings which do 
not have a distributed mass, are impacted in its middle point by the same 
body and have these characteristics: the light blue string is tangent to the 
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ideal string at its ends hence tension in these points is the same for 

both strings; the green string, which is not tangent to the ideal string in any 
point, has the same extremes and its middle point S 0 = 0 coincides during the 
movement and does not have also a distributed mass. Both delimiter strings 
shapes are always triangular during the loading and unloading stages, and to 
consider the contribution of the wire rope inertia the impact mass is redefined 
as a combination of the impact body mass, M c , and the wire rope mass, 


M 


WireRope 


Furthermore, it is assumed that the light blue string has small 


displacements while the green one has large displacements to be able to 
model the impact body movement. This does not mean the green string 
tension approximates the.yt^f^i0iR*%l>#ie impact point of the ideal string 


because its triangula 
point. This green strir 
body. Thus, the geom 
mathematical model 
attached to two linear 

■■■ Left spring tension level. 

1 1 Right spring tension level. 



ts deformation is larger at that 
the movement of the impact 
It blue string can be written as a 
the movement of a mass M 



due jv 


i T 


a 


a 

TUT T : I Z 

T 1 

*1 

1 

I 

■M 

’y 

^ R; * 

—►! 




Figure 15 Analogy between the delimiter model of tension at the ends of the string and 
a model of ideal springs. T, is the springs tension, a, is the springs slope angle, F, is the 
reaction force over the impact body, M, is the mass of the impact body, Y R , is its 
displacement, x, y are the coordinates and L, is the anchorage distance. 

Idealization of linear springs does not consider its real geometry and the 
model says the force felt by the ends R iy R 2 and mass M, T, is 
proportional to the displacement y M and a function of time. This is, 
transmission of force from the impact point and the string ends is 
instantaneous and furthermore the stress-strain relationship is always linear 
for the loading and unloading stage as shown in figure 16. 
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Figure 16: Stress (a)-strain (e) curve of a loading 
and unloading cycle of a perfectly elastic wire 
rope. 

The mathematical model of this movement with small displacements is: 



F= 


The last equation, (19), is valid only if tension T is uniform all along the 
wire rope or string. Now, wrCTrro^puint of reference in the impact over an 
ideal string the tension levels are as follows: 


T=AE 


Vj4+L 2 /4—L/2 1 


1. Force is a function of position and time: T(s 0 ,t) , with 
S 0 G [—L/2,L/2]. 

2. Because anchorages are considered perfectly rigid then on this points 
the tension is duplicated (reflection in a rigid frontier). 

And here is presented a geometrical model to show how forces change in 
a delimiter string which has, in every moment, a triangular shape and a 
concentrated mass on its middle point. Changes are presented for some 
multiples of the time to reflection r K =0.5L/C n with C n = ^Elp \ 
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1 W> 



► T(-L/2,4t R ) T (0,4g ^C„ T(0,4t R ) T(L/2, 4t R ) -* 

s ^ e tension level L- Overlap and add each other along the wire rope 
PvTI Right side tension level J 


Figure 17: Change in tension along the points of a wire rope idealized as 
a string impacted in its middle point and taking a triangular shape. M, is 
the mass of the impact body, a, angle of each right side of the string, F, 
reaction force over the impact body, T, distributed tension over the string, 
Cm, is the longitudinal stress wave velocity front velocity. The wire rope is 
idealized as a string without mass and with a “delimiter movement”. 


39 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Peru. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 




























































BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DE INFORMATICA Y COMUNICACION 


According to figure 17, it is not until time lr R that the string has 
stresses in all its points and, in this case, are symmetric with respect to the 
middle point. But it is not until time 2r R that tension can be considered 
uniform in both halves, this is, leveling of tension happens after a cycle of 
reflection when the wave front, with a zero stress level, is back at the impact 
point. Then, the same happens at time 4t R , next at time 6t R and so on. 
Thus, it can be say that from time t = 2t R and after time intervals At = 2t R , 
tension in the wire rope can be considered uniform along it. Hence, equation 
(19) is valid after the same time intervals 
2it R ; i = min ( t total /t R ), V i E IN — { 0 } . But as time to reflection t R is 



equivalent to travel LI2 ?,nd because the impact point coincides with the 

% 

after a distance 2C n (it R ). After 

W\ 


middle point then the ? 


this, the next expressi 


F=2AE 


= AE 


C u (*«*) 


Another way to express this models “the transverse movement of a 
material point attached to tWtP^iw§s as in figure 15 is valid only after time 


intervals At = 2t R ”, this is, each time this interval passes the next happens: 


r(-L/2,it R )«r(0,it R ) ; Vi G IN-{0} (21) 

As the time interval At = 2t R is very small compared with the time for 

loading and unloading, then, it can be considered that equation (21) is valid 
for these two stages. Hence as it R also comprise all this time, varying i is 
possible to approximate the movement in every instant by writing it R **t and 
the following: 


F=AE—^ = AE — 
C ii Mr C a t 


( 22 ) 


The physics expressed with geometry before can be translated to a 
natural language using terms from the method of changes used by Newton: 
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“Each change of the reaction (contact) force F between the wire rope 
and the impact body, this is, F+d t F(At) causes a change in the wire rope 
tension, which is T+d t T(At) . Then this value travels from the impact point 
through the string with the velocity of sound C n (for that wire rope 
construction and material), after a time it reachs and reflects at the ends, 
considered perfect rigidly here. There the wave duplicates its amplitud, 
reverses its sense of motion but not its nature (it is still tension 8 ). This means 
the waves from each side of the string, here equals because the impact is 
over the string middle point, after a time, overlap and add when they 
coincide along the string points. But, as the time taken by each change of T 
to cover all the string is a]sc«> then is considered that the level of 

tension after time inte v 
instant t. Also, bedu oar 


considered this results 


In equation (22), 



liform and by extension to every 
friction nor some damping is 
otion” 


is that point S 0 = 0 , or the impact 


point “see” changes trailing a distance^C^r and not simply LI 2 . While 

fj ' r Au ns 

division by 2 of L here is turnetJ into multiplication by 2 of the force which, 
geometrically, is the duplicative effect of reflection at the perfectly rigid ends, 
and it is also equivalent to travel half of the string back and forth by the wave 
front or the entire length in once as the impact is at the middle point. 

Now, it is possible to apply Newton’s laws to model the movement of the 
material point. And the next differential equation is obtained: 


dfy R y R 

M- lZ ^=-AE- r 


d.r 


Cut 


(23) 


But as was observed with the use of Yu. A. Dem’yanov’s method and 
the work of N. D. Cristescu, its equations can not be integrated from time 


8 Reflection of a wave is equivalent to the overlapping and addition of two identical waves, being of 
tension or compression, traveling to each other and when they met, the points reached by both are 
stretched or compressed by both waves and as they are equal then the amplitude is duplicated. 
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zero. The same happens here, thus,the movement is considered from time 
2 t R [66], [37], And equation (23) becomes: 


d]y R 

d.t 2 


AE y R 
M L + C u t 


(24) 


This equation, (24), is of the Emden-Fowler type and according to K. S. 
Govinder and P. G. L. Leach, this type can be integrated and has an 
standardwhich is [68] 



T^=f(t)y" 


In equation (24), function f and coefficient n are: 


Then, knowing tf 
the lineal solution of ©goal 
following shape: 


(26) 

i approximate solution analog to 
established, and should has the 


y R .a P ro X B 6 t) 

Its initial conditions are: 


; d ^f( t =°)= v « ( 28 ) 

The mass to consider is M=M c +M cable . And an approximate solution 
for the displacement (model of small displacements) of the light blue string’s 
middle point is: 




aprox . 




EA t 
M V L+C n t 


(29) 


Equation (29) describes a resonant system which amplitude grows and 
its frequency decrease continuously with time, this is, only represents a 
solution for a first loading stage. And, to test if this is an approximate solution 
it is replaced in equation (24) as follows: 
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d 2 t y R . 


aprox. 


u 2 y R . 


aprox. 


-u\y R . 


aprox 




4 L + C n t 


L+C n t 


Oh 


t/ 


4(L + C w t) 


9{t) 


U 1= ^M ; U 2 =lJ** 

1 2 ] EA 2 ] M 


( 30 ) 


Where: 


j(0= 


1 

c 1 

11 + 

L + C n t 

4 U\ 

4 AE 



C„t . 

C 2 

'-‘ii 

t 2 

2 

4 

L+C u t 1 


M(L + C„t) 4lfeC„t) 


: COS 


U 2 t 



L + C n t 


(31) 


(32) 


(33) 


(34) 


If every constant aie evaluaieu j(t) and g(t) for r<l then: 
L^iMlAp , as E has an order or 10 u [Pa] , M and order of 10 3 [Kg] , A 
of 10“ 4 [m 2 ] , p of 10 3 \Kglm 3 ], C n of 10 3 and V up to and order of 10 1 . 
And with all these: 


j{ 0>i ; 


4 AE 


g(t)**0 


(35) 


M(L + C u t) 

And j(t) , with the experimental dimensions and suppositions its value 


tend to be 1.1091. Hence a correction factor can be established for y R aprox 


9. 


J 


(36) 


/ Corr. 


— lim 


= 1.1091 


9 This calculation can be probed using QtiPlot and the attached fdes. 
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Thus, the corrected displacement of the string’s impacted point with a 
delimiter movement and small displacement is: 


, . dfy R 

y R.corr. rCorr. y R.aprox. S 

d t t 


-u; 


L + C n t 


&Corr.y R. 


aprox. 


(37) 


y R.Corr. TCorr. 


= <Pcorr.^y-^-^JL+C n tsenl2^ 


EA 


t 


(38) 


2 V EA 11 | V M ylL+C n t\ 

Then the equation that models the force exerted at the ends or 
anchorages during the loading stage or deceleration of the impact body is: 


T C =AE 


4> 


Corr. 


This equation is 
reached or only for the 



+ 1-1 


(39) 


tact up to a maximum tension is 

tr 


Hence, a model for f he axial^stress (as an average value for all the 
wires) at the ends is: 


o l =E 


00 


'Corr. 


\f A ^ 


tsen 


1 


2 

■J§' 

(-f 
v 2 ’ 

+ 1-1 

V L+C n t 


(40) 


for 0 ^tKt c mAx i 

The last two equations represent a resonant system and only model the 
loading stage. This is why time is restricted up to the maximum tension is 
reached or the mass M is completely stopped. Then, just an instant later, 
the unloading stage begins. Thus, to model this process other equations are 
needed. To take into account the unloading stage, the elastic modulus of a 
wire rope has to be analyzed again. This parameter changes with every cycle 
and level of load or during unloading, and in general, in this stage its value is 
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bigger than when loading. Due to this dependency and knowing that a way to 
have a predictable design is having constant properties, there is a reason 
why it is recommended to load a new wire rope ten times before measuring 
its modulus of elasticity. But all these properties are given for quasi-static 
applications with little ratios of deformation e:10~ 5 tol0~ 1 mm/mm/s . And this 
is not enough to even classify a movement at the first scale of the 
Characterization of Weihrauch ( e with an order of 10° mm/mm/s ) [43], [69], 
But high ratios of deformation are present during the transverse impact of 
wire ropes as the ones in the references used to validate the models here. 
There, the ratios can be up to e:2mm/mm/s . And this, of course, affects 
how the elastic modulus charpg^y^ji a wire rope is loaded. Observing the 
references results (us 1 I jn be say that as with a simple 


carbon steel bar, a I 
changes the elastic nr 
almost static load. The sarr 
structure would be me 



on makes it more resistant and 
si than when it is loaded with an 
r a wire rope and in this case its 
ading of a wire rope is described 

a new one 
being loaded and unloaded (and 


using the curve in figure 16, but ever^Sfre rope and specially 
presents a permanent deformation after being loaded and uni 


presents a permanent 
more if there is plastic deformation). This is why a better approximation is 
given in figure 18, where a wire rope is unloaded from a maximum 
deformation e Mdx to a final e f >0 with a higher elastic modulus, as it is near 
to be constant and more stable. 



Figure 18: Stress (a)-strain (e) curves of a 
loading and unloading cycle of an idealized 
wire rope. £ M ax, £f, are the maximum and 
minimum deformation respectively. 
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Also here is supposed that the unloading stage begins after the impact 
body reaches a null velocity and all its momentum has been transferred to 
the wire rope which, at the same time, has a maximum tension or a 
maximum deceleration. Then after this instantaneous tension the inverse 
begins to happen, this is, all the stress accumulated by the stress waves 
reflection begins to decrease in the same way as it increased but with more 
speed because the elasticity modulus is bigger. With geometry, the latter is 
expressed in figure 19, where an space without gravity or resistance medium 
is represented along with two perfectly rigid strings with uniform and identical 
properties are under longitudinal impact by a punctual mass which has an 



Figure 19: Illustration to demonstrate the symmetry with respect to time of a 
loading and an unloading process of two identical perfectly elastic strings and 
with the same maximum and minimum tension levels T. Both are longitudinally 
impacted by a punctual mass M, i, represents equally spaced time instants, 
t n , represents the period from an minimum to a maximum and back to a 
minimum, g, is gravity, x,y, are the spacial coordinates, S, are the curvilinear 
coordinates. 
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Figure 19 can have this interpretation: first, unloading of string A and 
loading of string B is expressed in the same time interval re[0;t n ], where 
string A or B reaches its maximum tension at time t n (the same for every 
point S 0A ) when loading or its minimum after unloading, this is, in this case 
t c.Max ' s how longer the loading stage last, and also coincides with the time 
that the unloading stage endures: t D Mdx . Now, if this time interval is divided 
into equal parts At = t n /n , then At=t i -t i _ 1 ; l</<n,Vi,nelN-{0} , thus 
it is possible to write the unloading tension of string A as a function of the 
tension in string B while it is loading for each time r, and every of its points 



(41) 


iAt and as i'Ar = r ; then 


t n _.=t n -ti , thus tension of string A while unloading can be expressed with 
T A {ti)=T B {t n -ti) and because the process is continuous, taking the limit, 
results in: 


T A (t) = T B (t„-t)=T B (-(t-t„)) 


(42) 


This is, the loading and unloading process are symmetrical functions of 
time with the values of tension T B transferred a time r n thus both stages 
take place at the same time. This is also compatible with the symmetry 
principle of time of Emy Noether, which says (basically) “ what is before is 
after” and this way conservation of energy is established which in the case of 
strings under loading or loading is conservation of elastic potential energy 
[70], Then, loading tension as a function of time and its symmetrical process 
of unloading of string A and B can be expressed in figure 20 as follows. 
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T D (t)=T csim (-(t-g) T D (t)=T asim (-(t-(t DMa +t CiM J)) 



^D.Max"^C.Max t 

Figure 20: Repress^"' nsion Vs Time of a loading and 

unloading process f>W m i. according with the wire rope 
behavior of figure • iviJNf j rapft . i'm tension while it is loaded or 

the beginning of the r^^e^^r^this is the maximum tension at 

the beginning of the n the unloading finish, T D , level 

of tension while un Tc.sim., level of tensions while 

unloading as a symmc' r ic function of the loading stage. 

Now, the unloading tension T D (t) of any ideal elastic string and with 

constant properties can be expressed as a symmetrical function of a loading 
process of the same string: 


^d( 0 — ^C.Sim. ( (f ^C. Max) ) (43) 

This equation is valid at the interval t c Mdx <t<t c Mdx +t D Mdx being the 

instant t c M6x , where the loading stage ends and the string tension reaches 

its maximum. But here (process of figure (20)), as both strings are the same 

then t c _ Mdx =t D Mdx . The same figure shows how a symmetric function can be 

translated to any instant, hence, both stages can happen simultaneously by 

using expression -( t-t C Mdx ) and the unloading stage can happen just after 

the other is loaded using -(t-2t c _ Mdx ). 


48 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Peru. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 











BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DE INFORMATICA Y COMUNICACION 


And because strings A and B have no mass or bending resistance, 
then, its state, in any moment, is equivalent to that of a string with a “delimiter 
movement”, with a triangular shape, small displacements and tension equal 
to the tension at the ends of a transversely impacted ideal string like in figure 
14. And this way is possible to continue approximating the transverse impact 
over wire ropes. 


But, as the change of modulus of elasticity after the loading stage is 
discontinuous this does not imply necessarily that changes in tension are 
also discontinuous. The maximum tension reached after loading 
T c . m =MAX{T c ) with T c from equation (39) has to be equal to the 
maximum tension reachedy&frer a unloading process symmetric to one of 
loading in the same w ggfrrtpf ^ing modulus, this is: 

T c 

In a real applicat 
possible to say that 



Max) (44) 

ulus is progressive but here its is 

ading modulus of elasticity are 


related with E D =R 78 E where J? 78 >1 ,vtd a constant. And to show that 




0 


equation (44) is true, equationused for both stages: 





As E d —R 78 E , equation (46) can be expressed as follows: 


T D.Mdx.— AE 


R 


i y R . 


Corr AC. Max) 


781 




R 


78 \ 


LI 2 


+ 1-R 


78 


COR y R _ Co rr. S( E D. Mdx)~ y R.Corr.(C. Max) 


/V R 


78 


(47) 
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In equation (47), y R Corr S(£D) , is the displacement during the unloading 
expressed as a symmetric function of the loading of the same string with its 
unloading modulus of elasticity E D and with the model of small 
displacements of equation (38). And because y R Corr (c Mdx] , the maximum 
displacement reached during the loading stage, is a constant then 
(y R .corr. { c.Ma X )/ L l 2 ) 2 =R79 is a constant too, then: 


V ^79"*" 1 1 — V^78 V^79+^R78 i?78 (48) 

If roots in this relation are expanded using the first terms of a McClaurin 
series then the maximum tension reached is the same with which the 


unloading stage begins: 


Knowing this, it 
tension at the ends d 
the loading of the sam opc 


T D — AE D 





M 


W=R 7 o ( 49 ) 

5 the model to approximate the 
ige as a symmetrical function of 
unloading modulus E D : 


sen 



+ 1-1 


(50) 


Valid in the time interval t c . MAx <t<t D . MAx +t c . MAx . Being, t DESC : 


tDESC~{ D.MAX + tc.MAX ))) ’ tc .MAX < ^^^D.MAX + ^C .MAX ( 51 ) 

In (50) the symmetrical functions expressed as functions of (—t) have 


to be transferred to t 


C. MAX 


because it the unloading stage begins at that 


instant, and this is why t c MAx +t D mAx 10 is added to -t. 


10 In practice, determination of these times can be done with easiness using some software. 
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Before writing a model for the movement of the impact body, the model 
of the displacement during the unloading stage Corr D °f the delimiter 
string with small displacements has to be established. According to equation 
(50) displacement during the unloading stage as a symmetrical function of a 
loading process with module E D is: 


R.Corr.SlED, 


= <P 


00 


Corr. 


I M 

1 \e d a 


V L + C n t 


DESC 


sen 


l 


E n A 

211 


DESC 


L + C u t 


DESC 


(52) 


And as E D is bigger than the one during loading then the maximum or 
initial value calculated with (52) is smaller to the final or maximum 


displacement reached d' 
this is compatible wit|| 
in figure 18. Although, 
delimiter string middle 
during the unloading t 



9®' y R.Corr.SiED.Mdx^y R.Corr .(C. Max) ’ 

in the unloading stage showed 
)t model the displacement of the 
how much stretch has the string 
it the actual displacement (small 


displacements model) in this stage in aigay that the total displacement is a 
continuous function, the changes of coordinate of the middle point caused by 
its velocity has to be added, this is, the derivative of equation (52), 
y ' R Corr S(ED) has to be integrated and added: 


y R. Corr. D~ y R. Corr.(C.Mdx) + j 


(53) 


V D , is the velocity during the unloading stage, it goes from zero to a 
maximum when the string tension is zero again. This velocity is: 


V n = 


D d t t 


y R . 


Corr.S(ED) 


(54) 


Where: 


R. Corr. S[ ED 


=y 


R.CorrAED 


HHtD. 


MAX + tc.MAX 


))) 


(55) 
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And the displacement (small displacements model) during the 
unloading stage is: 


t 


t 


J* Vod t t ^ ^ J y R.Corr .S(ED)dtt y R. Corr. S(ED)(0 T R.Corr. S(ED){^C .Max) (56) 

tc.Max ^C.Max 

If time t c Mdx is replaced into equation (55) then the term after sign 
minus in equation (56) is equal to the maximum reached during a symmetric 
loading process with an elasticity modulus equal to E D , as follows: 



y R. Corr. S(ED){^C. Mdx )~y R. Corr. (ED){^D. Mdx ) — T R.Corr .(ED). MAX (57) 

And replacing this result into equation (53) then the displacement during 
the unloading is: 


.Corr .D " =y R .Corr .(C 

Valid in the ti 


(* (* d.mAx + ^c.mAx 


))) (58) 


This way the model of smlfll ^replacements is continuous in every 
instant no matter if the elastic mod^^ehanges discontinuously. 


In sum, equations (38) and (58) model the small displacements of the 
middle point if an string which movement delimits that of an ideal string 
impacted on its middle point and has the same tension than the ends of this 
ideal string. 


4.2 . Mathematical model of the impact body. 

Now, it is possible to model the kinematics of a body idealized as a 
material point that impacts an ideal string. This begins with the delimiter 
movement of a string with large displacements represented with the green 
string with triangular shape in figure 14. But first, as the lineal case has the 
following relation: 
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F= 2 


AEy, 


y r _ 

con —— » 


L/2 ’ —' LI 2" +/, <60) 

This means the angle a is assumed small and this way the stretch in 
the middle point S 0 (L/2,t ) can be approximated with the value y R . But here 

a better approximation is VL 2 / 4+y 2 M -Ll2 . It is important to emphasize that 
the tension of the ideal string is not that of the green string because the 
difference in shape and stretch. Then, and alternative to find a mathematical 
model is to establish that the stretch of the string with large displacements 
can be approximated for the loading stage y M c , and unloading stage y M D , 
as a function of the displacements of the string with small displacements 
y R Corr . Thus, for both sta<js&: 

V L 2 /4v 
With this, accek 

mass M while the str 



% "2 (61) 
2,t) or of the impact body with 


^ Tm.C _ „ AE I y R .Corr. Tfi.Corr. ' 

i/ = 4* ’ter, 


Vl 2 /4 +> -i 

for 0<t<t c mAx 

And the acceleration while unloading is: 


(62) 


y m .d _ D 


df 


M 


y R . 


Corr. D 


y R . 


Corr. D 


I 


Ll 2 Jl 2 /4+v 2 

V L ' y R .Corr.D i 

f° r t C .MAX<t^tD. M AX + t C .MAX 
Velocity of the impact body M in the string loading is: 


(63) 


d t t ~ V0 ° J 


f 2 AE 

y R.Corr. 

y R. Corr. I 

J “ M 

0 

L/2 

^l 2 /4+ y R Corr J 


d t t 


for 0<t<t 


C.MAX 


(64) 


53 


Es/a obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Peru. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/Hcences/by-nc-sa/2.5/pe/ 




























BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DE INFORMATICA Y COMUNICACION 


Velocity while unloading is: 


d t y m.d. _ 

d t t 


f 2 AE 

y R.Corr.D 3^ R.Corr.D | 

J 1 M 

mAx 

112 <jL 2 /4+y 2 R Corr D J 


(65) 


Para el intervalo t c t^t D max* t Cm max 
In the equation above, V F c , is the velocity at the of the loading stage 

and is not necessarily zero in this model, because is approximated, this 
is,when the velocity is effectively zero, the calculated force is not necessarily 
maximum at the anchorages. V F c can be calculated with: 



Displacement 


string loading is: 


y M.C. — V oot 


y r. 


Corr. 


L 2 /4+y 2 R 


. Corr. i 


for 0 <t<t CMAx 
Displacement while the string unloads is: 


( 66 ) 


(67) 


yM.D-Vp.C 


-J 


L C.MAX 


J 


AE, 


M 


y R . 


Corr. D 


y r. 


Corr. D 


Ll 2 Jl 2 /4+v 2 

V L J R.Corr.D 


I d t t 


L C .MAX 


d,t ( 68 ) 


f 


f 0r MAX + t C . MAX 

And finally the force of reaction F between the wire rope and the 
impact body is: 


F C =-2AE 


y R. Corr. 

i 

y R. Corr . 

L/2 

'J L 2 /4 + y R Corr 


for 0<f <r c M ^ x 


(69) 
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And while unloading is: 


F d =-2AE d 


y R . 


Corr . D 


y R . 


Corr . D 


L/2 


V L 2 /4+ 


( 70 ) 


' yR.Corr.D J 

Para el intervalo t c M ^x<t^t D M ^x+t c 
Equations from (62) to (70) can be integrated numerically with a 

program like GNU Octave or using QtiPlot. 

Nd as a commentary, even with the discontinuity presented by the 
change of the elasticity modulus, all constructed here can be used directly as 
frontier conditions in a numerical method to solve the partial equations (8) 
and (9) with the only restriction that the string have to be elastic. This method 
can be similar to that .'^ffeje^^Nf^^^tescu for a string with a free end 


but now developing th.~ 


4.3 . Mathema 
the impa 
rope. 



he author of this work [37], 

he internal stresses in 
sversely impacted wire 


The internal stresses arise-^^^^the contact between the wire rope 
and the impact body. There emerges a force of reaction F , which is 
transferred to all the wires causing not only a resultant longitudinal tension 
but shearing, bending and torsion stress and also contact stress between 
them. This way the wire rope presents higher stresses around the contact 
zone compared with what the anchorages support. Thus, the reason of this 
modeling is to include these additional stresses in the wires and then 
combine them with the longitudinal stresses of the wire rope and end up with 
a criterion to determine a security factor for the wire rope analyzed. 

To model these stresses there are some options, as: 

1. It is possible to consider the complete wire rope geometry and that of 
the impact body: the curvature of the surfaces in contact, of the wires, 
the force between each wire and other interior components. But this 
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leads to a very elaborated geometry and represents another research 
work. 


2. Another alternative is to use only the surfaces in contact but this way 
the wires and strands re-arrangement is not taken in count. 

Then, a more viable option is to work as G. Rebel and R. Verreet did to 
study the surface wear of wire ropes under pressure when are rolled over 
drums. They modeled the contact stresses using Hertz theory, taking the 
nominal diameter of the wire ropes, a Poisson coefficient equal to 0.3 (as for 
steel) and the wire rope elastic modulus (smaller than that of steel) [74], This 
way it is possible to approximate the stresses along the z axis(figures 21 
and 22), and only thes° f he maximums around the contact 

zone and represent 
impacted 11 . Two situal 


in contact with the eq 


D 2 , their axis are at 



' that arise when the wire rope is 
A steel cylinder with diameter D x 


d of the wire rope with diameter 
:he second case is a rigid plane 


front body of length n ;n contact with^e wire rope equivalent cylindrical 
solid with diameter D. These cases are in concordance to what were used 
in the experiments and simulations used here to validate the models. The 
equivalent diameter of a wire rope is found using its metallic area: 


£> 2 =£>=2Va/jt (71) 

The case of two perpendiculars cylinders is done assuming that the 
contact area is circular, next, an equivalent radio for the two cylinders in 
contact is found as follows [71], [72], [73]: 


(72) 

Where H 1 and H 2 are calculated using the diameters of both cylinders 
in contact D, y D 2 : 


11 All the stresses in the contact zone and around it can be determined using the procedure and solving 
the elliptical integrals presented in references [71], [72], [73], 
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H 1+ H 2 = ±[2/D 1+ 2/D 2 ] 


( 73 ) 


H 2 -H 1 =|[4/Di+2/D2-8/(D 1 D 2 )] 


1/2 


With (74) is calculated the contact area’s radio: 


(74) 


a= nWU-vjAuJ 


k 1/3 


(75) 


In (75), (i is the shear modulus ^=£72(1+ v). And the stresses along 
axis z are: 



0 XX Gyy 


p 0 , is the 


i-i'] 


(77) 

over the circular contact area 


and depends on the maximum reaction force F mAx 


Po— 3F MAx/2-Jta (78) 

The equivalent perpendicular cylinders are as shown in figure 21: 



y\ 


Figure 21 : Cylinders of diameters Di, and D 2 in contact with its 
axis at 90°. x,y,z, are the coordinated axis with origin at the 
center of the contact zone. 
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For a cylinder of diameter D in contact with a plane rigid body the 
contact is a rectangular area formed over the cylinder with a width b : 


rj 1/2 

b = (^FD[(l-v)/E]) 


(79) 


The stresses along z axis in this case are: 



(80) 


a yy =-2p 0 v((l +z 1 lb 2 ) lll -zlb) 


(81) 



(82) 


Here, p 0 , is also the Tl ^;SQi£/j T /i|jni pressure reached over the contact 



Figure 22: A cylinder with radio R and a rigid plane of length n in 
contact, x, y, z are the coordinated axis with origin in the middle of the 
contact zone. 

Now it is possible to establish a model for the stress state around the 
contact zone during the transverse impact. This is represented in figure (23), 
there is included only the compression stresses which arise by the contact 
along z axis, and are also principal stresses <y xx ,o yy ,o ZI calculated for the 
maximum reaction force F Mdx =Ma Mdx , there is also added the stresses 
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caused by the tension along the wire rope or the y axis, here 
cr L.MAx( s o = 0,t) which can be calculated or measured. 




a yy +a( s o=°, t ) 


Figure 23: Maximum Rfiitfe^lfcf/t^ess state along axis z inside a 
transversely imp At-'’ . io th Q ma^tinn f<~»™ 

between the imp 

I 


maximum tensiti? 
o xx ,is the princi| 
parallel to y axis : 
longitudinal stres 



e rope surface, T m a X . is the 
s rope at its anchorages, 
x axis, a yy is the principal 
ress parallel to z axis, a , 


The foregoing le imp-® erion f° r the impacted wire rope, 

that could be establishectaising the Vom^Tses criterion (by analogy between 
a ductile material and a wire' ropJe which structure rearrange while being 
loaded) to relate the internal stresses by contact and the longitudinal 
stresses, and then to compare results with the material strength limit and to 
determine a security factor. It has to be noted that the longitudinal stress at 
the point of impact (or zone) is not directly determined here but as the 
maximum value is needed, it is approximated using the maximum reached at 
the ends. Thus, using the Von Mises criterion along the z axis with origin at 
the impact point, the stresses are: 


oV\ j ^ ^ Tv O'l.max) (°yy max ( rr 7/ r AV ) 


(84) 


And this can be compared with the maximum strength of the wire rope 
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5 . Validation of the mathematical models. 

The experimental data used to validate the model of longitudinal forces 
and the movement of the impact body were collected in the experiments of C. 
Stolle and the simulated data were obtained in simulations done by Tytko et 
al y P. Gospodarczyk [31], [29], [33], 


5.1 . Validation of the model of longitudinal stresses 
and the model of the impact body’s movement 
using the experimental data of C. Stolle. 



In this part, models are applied using properties and conditions as used 
by C. Stolle on its experiments. Tension at the ends of the wire rope while 
this is loaded or unloaded are ca!cu..,ted using equations (39) and (50), 
acceleration of the irr^ v d with (62) and (63), its velocity 

with (64), (65), (66), a 


C. Stolle did 6 


impact on its middle 


equations (67) and (68). 

el wire ropes under transverse 
xperiments n° 1, 2, and 4 12 are 
taken. Wire rope characteristics necessary for mathematical calculations are 
in table 3. 


Wire rope type 

Number of strands 

Total number of 
wires 

Nominal 

diameter 

3x7 

3 

21 

19 mm 

Wire rope 
metallic area 

Weight per unit 
length 

Steel density 

Breaking 

strength 

1.54508E-04 m 2 

1.22803E+00 Kg/m 

7.80000E+03 Kg/m 3 

1.80100E+05N 

Elastic limit 

Elastic modulus in 
an static test 

Elastic modulus of a 
pre-tensioned wire 
rope 

C. Stolle 
recomended 
elastic modulus 

1.10400E+05N 

7.99500E+10 Pa 

1.16000E+11 Pa 

1.04000E+11 Pa 


Table 3: Physical characteristics of a wire rope type 3x7 used by C. Stolle in 
his experiments [31], 


The steel wire rope anchorages were conical sockets filled with epoxy 
resin and joined to a steel base which let the rotation of the anchorage 


12 C. Stolle, mentions that in his experiments 3,5 and 6 occur a lost of data and technical problems 
and the data were not appropriate for his simulation and are not for this work too. 
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around a vertical axis. This base, in turn, had a reinforced concrete base as 
shown in figure 24. 




0 19 mm 


Figure 24: Photograph r 
socket used as ar.:!. 
wire rope used by 
experiments [31]. 

C. Stolle, alr~ 

frontal part inclL::' _ 

a way that it was i. 

movement and no l 


0 3.061 mm 


lure 25: Cross section of the 3x7 
'i rope used by C. Stolle in his 
iments. 

for his experiments, in its 
'7..imandA36 (figure 51) in 

.i the wire rope during the 

_; with the car (impact body) 

movement. A rail was use6Md<#>EStt^^the tires and the car straight and 
perpendicular to the wire rope while it was impulsed by a pickup up to the 
wanted velocity of impact. Figure 26 shows the car contacting the wire rope: 



Figure 26: Impact position of the car and the 3x7 wire rope 
on the transverse impacting experiments of C. Stolle [31]. 
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In this section the developed models are used and its results are 
compared with the experimental data of C. Stolle, and as him also did a 
numerical simulation this results will be used to compare the models. In his 
experiments the wire rope length was larger than the distance between 
anchorages: 46.2 mm more in experiments 1 and 2 and 50.8 mm in 
experiment 4. The wire rope was brand new and its shape before the impact 
had waves due to its winding over a drum and gravity. C. Stolle did included 
the additional length and the non-straight shape in its simulation and also 
added factors to model the internal friction, bending resistance, torsion and 
plastic deformation. 


But, in the other hand^t^hatbematical models here do not take into 

to consider into the models are 


account those factors # 
only the properties gft> 

V 00 , its mass M c , tf 
figure 27, and tables £ 



velocity at the instant of impact 
nchorages L as represented in 


5- f l 772.0; 772; 769.3] [ Kg] 
SaaSPC- [6.625; 6.75; 10.63] [m/s] 


Rigid anchorages 


Wire rope 3x7. 



Figure 27: Idealization of experiments 1,2 and 4 of C. Stole. L, is the 
distance between anchorages or ideal wire rope length, M c , is the 
impact body mass, Voo, is the velocity of impact, x,y are the coordinated 
axis and O, is the middle point of the wire rope and coordinate origin. 
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5.1.1 . Modeling experiment 1 of C. Stolle. 

Value of entry variables and the parameters considered here are in 
tables 4,5 and 6. Results are presented in figures 28, 29, 30 and 31 and 
tables 7,8,9 and 10. 


Impact velocity 

Mass of the car 

Distance 

between 

anchorages 

Initial wire rope 
tension (New wire rope 
without pre-tension) 

6.62500E+00 m/s 

7.72000E+02 Kg 

8.61100E+00 m 

0.00000E+00 N 

Table 4: Entry variables to model experiment 1 of C. 

Stolle [31], 


Impact body 
(car) mass M c 

7.72000E+02 Kg 




Wire rope 
metallic area 
1.54508E-04 m 2 

M= M + M hl 

c cable. 

7.82575E+02 Kg 
Table 5: Parameters to 
loading stage. 



Meal cable length L. Elastic ^ odulus 


%• OOE+OOm 


1.04000E+11 Pa 


d velocity in Wire rope mass 


steel C 


M 


cable 


48E+03 m/s 1.05746E+01 Kg 

^2 ^c. MAX 

89E+02 m 05 /s 1.84057E-01 Seg 

lent 1 of C. Stolle for the wire rope 


Impact body 
(car) mass M c 

Impact body 
velocity V 00 

Ideal cable length L. 

Elastic modulus 

E 

7.72000E+02 Kg 

6.62500E+00 m/s 

8.61100E+00m 

1.68411E+11 Pa 

Wire rope 
metallic area 

Wire rope steel 
density 

Sound velocity in 
steel C M 

Wire rope mass 

^cable 

1.54508E-04 m 2 

7.80000E+03 Kg/m 3 

4.64662E+03 m/s 

1.05746E+01 Kg 

M=M C+ M cab|e 

u, 

u 2 

^C.MAX 

7.82575E+02 Kg 

1.81660E-02 m 

3.64693E+02 m° 5 /s 

1.84057E-01 Seg 


^D.MAX 

^C.MAX + ^D.MAX 

1.44601E-01 Seg 

3.28658E-01 Seg 


Table 6: Parameters to model the experiment 1 of C. Stolle for the wire rope 
unloading stage. 
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Table 7: Comparison between the maximums of the tension at the anchorages of 
the transversely impacted wire rope of the Exp. 1 of C. Stolle. 
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3.28090E-01 Seg. 3.50000E-01 Seg. 3.00000E-01 Seg. 

Table 9: Maximums return velocities of the impact body, calculated 
with the measured and simulated in the Exp. 1 of C. Stolle. 
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1.84095E-01 Seg. 1.71030E-01 Seg. 1.93400E-01 Seg. 

Table 10: Maximums accelerations of the impact body, calculated 
with the measured and simulated in the Exp. 1 of C. Stolle. 
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5.1.1.1 .Discussion of the modeling of experiment 1 
of C. Stolle. 


Before using the developed equations, all the parameters have to be 
determined. Values of the material point mass M , its impact velocity V 00 
are not a problem, but as the wire rope has slackness between the 
anchorages, this requires a decision to make about what is the value of L , 
and here is taken as the anchorages distance knowing its difference with the 
real length is less than 50 mm. Another decision to make, and more 
important, is what value can be chosen for the elastic modulus knowing it 
changes with time and load; about this, C. Stolle recommends a value of 104 



Gpa for a new wire rope and 115 Gpa if it has a pretension but only to 


calculate the sound v .• s >' ^ 

/fy 

3x7 wire rope has a rr 


not give precisely wh; 
stage. Thus, in this pc 
to use in the develope 


"■%js a mean value because a new 
% its first load test. C. Stolle does 
imulate the loading or unloading 
as to be made about what value 
g the aforementioned, a value of 


104 Gpa is used as a mean for the loading stage because the experimented 
wire rope is brand new. 


For the unloading stage the next reasoning is followed: finishing the 
loading stage its structure is more or less under a pretension, more 
stabilized, therefore its modulus would be approximately 116 Gpa, as it is 
recommended by C. Stolle, but, with the high deformation rate influence and 
as the elastic modulus in the unloading stage has to be bigger than in the 
loading stage then its modulus will change from 79.9 Gpa to 116 Gpa at the 
end of the loading stage and will change in the same proportion after 
beginning the unloading stage, this is, from 116 Gpa to 168.46 Gpa as is in 
tables 6 and 13. 


About the results, in figure 28, the mathematical models predict a higher 
tension during the loading up to after the referential elastic limit at -0.12 S 
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(being 110 KN, for an static test). This happens because the chosen modulus 
is constant and bigger than a new wire rope has at the beginning and also 
influence the wire rope pull through the socket’s resin which also absorbed 
energy and caused that a minor tension is needed to stretch it at the 
beginning as the experimental data shows, followed with a fast rise of tension 
after 0.05 S. Then its structure stabilize more and its tension reaches and 
surpass the elastic limit while its elastic modulus continues to increase, but, 
then the wire rope begins to deform plastically and decrease until the 
unloading begins. And is in these two parts where the approximation is better 
because the wire rope structure is even more stable but differs again, as the 
predicted tension is higher after the elastic limit and continues this way up to 
its tension becomes zi r^pmis ' a* V. ause, while unloading, the elastic 


modulus is not constat 


Then, figure 31 
marked discontinuity 



'.own in figure 10. 


;ulated acceleration presents a 
jIus changes abruptly from 104 
nloading. Before this jump the 


Gpa, while loading, ^ 

analytical curve can be considered ae&a mean value of the acceleration 
during the loading stage. While for the unloading stage this curve presents a 
better approximation, as also happens with tension, and allows predicting the 
maximum acceleration reached. This is an important value to consider when 
designing vial barriers because it is a reference to know the forces exerted 
over the passengers in a car when they have an accident. 


The predicted velocity in figure 30, approximates with more precision, 
both, the loading and unloading stages and even the final or return velocity. 
This happens because, even if acceleration has a large discontinuity, the 
experimental acceleration is not smooth, has peaks for some hundredths of a 
second that after integration results in a smooth velocity curve. And as the 
calculated acceleration is as a mean value during the loading stage and 
approximates better while unloading, this results in velocity values very close 
to the real ones. The predicted zero velocity does not coincides with the 
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maximum force, as was pointed out before; this is a result of the difference 
between tension, in value and time, between the ends and at the impact 
point. The same explanation has the displacement prediction showed in 
figure 29, which is almost the same in every instant except at the end. 


In general the analytical predictions are not as much close to the 
experimental data as are the C. Stolle simulations in every curve. But some 
disadvantages to use his work are that he does not provide all the necessary 
data to use his simulation, and LS-DYNA, is expensive and not libre 
software. However, the developed models here can be used if the elastic 
modulus is close to a constant, or in general if the structure properties are 
uniform or constants, this ie.^fvfc^thematics is to be used then it is possible to 

stabilize the wire rope ; This could decrease its capacity 

^ 'is* 

to absorb energy bi# nas to have energy dissipators 


instead of allowing tl 
breaking point or inste 



} a plastic deformation near its 
engineered structure 
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5.2 . Modeling experiment 2 of C. Stolle. 

Entry variables and parameters necessary to model this experiment are 
in tables 11, 12 and 13. Calculation results are presented in figures 32, 33, 
34 and 35 and in tables 14, 15, 16 and 17. 


Impact velocity 

Distance 

Mass of the car between 

anchorages 

Initial wire rope 
tension (New wire rope 
without pre-tension) 

6.75000E+00 m/s 

7.72000E+02 Kg 8.61100E+00m 

0.00000E+00 N 

Table 11 : Entry variables to model experiment 2 of C. 

Stolle [31]. 

Impact body 
(car) mass M c 

7.72000E+02 Kg 

Impact body velocity ideal cable length m . .. . , 

w . Elastic modulus E 

V 00. 

6.7500QpS$®%' ,R.61100E+00 m 1.04000E+11 Pa 

Wire rope metallic 
area 

1.54508E-04 m 2 

v 0 > * 1d ve|ocit y in Wire ro P e mass 

steelc " Mcab,e 

7.8 48E+03 m/s 1.05746E+01 Kg 


M=M + M .. 

c cable. 

7.82575E+02 Kg 


u 2 

"89E+02 m°- 5 /s 


Table 12: Parameters 
loading stage. 


VMAX 

1.83612E-01 Seg 


nt 2 of C. Stolle for the wire rope 

w 


Impact body 
(car) mass M c 

-- 

Impact body velocity Ideal cable length 

V 00 L. 

Elastic modulus E 

7.72000E+02 Kg 

6.75000E+00 m/s 8.61100E+00m 

1.68410E+11 Pa 

Wire rope metallic 
area 

Wire rope steel Sound velocity in 

density steel C M 

Wire rope mass 

M hl 

cable 

1.54508E-04 m 2 

7.80000E+03 Kg/m 3 4.64661 E+03 m/s 

1.05746E+01 Kg 

M= M + M H 

c cable. 

u u 

W 1 2 

^c. MAX 

7.82575E+02 Kg 

1.85088E-02 m 3.64692E+02 m° 5 /s 

1.83612E-01 Seg 

^D.MAX 

t +t 

l c.MAX l D.MAX 


1.43980E-01 Seg 

3.27592E-01 Seg 



Table 13: Parameters to model the experiment 2 of C. Stolle for the wire rope 
loading stage. 
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Table 14: Comparison between the maximums of the tension at the anchorages of the 
transversely impacted wire rope of the Exp. 2 of C. Stolle. 
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Table 15: Comparison of maximums displacements of the impact 
body, calculated with the measured and simulated in the Exp. 2 of 
C. Stolle. 
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3.27590E-01 Seg. 3.00000E-01 Seg. 2.81300E-01 Seg. 

Table 16: Maximums return velocities of the impact body, calculated 
with the measured and simulated in the Exp. 2 of C. Stolle. 
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1.94500E-01 Seg. 1.88190E-01 Seg. 1.84094E-01 Seg. 

Table 17: Maximums accelerations of the impact body, calculated 
with the measured and simulated in the Exp. 2 of C. Stolle. 
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5.2.1.1 .Discussion of the modeling of experiment 2 
of C. Stolle. 

Here are used the same elastic modulus as for experiment 1 because 
the same new wire rope was used. Only the impact velocity changes a little 
percentage, from 6.625 m/s in exp.1 to 6.75 m/s, but this is enough to cause 
a larger difference between the calculated and measured tensions even after 
the elastic limit of 110 KN is reached while loading. But a better 
approximation of the maximum tension is also reached: in experiment 1 the 
model predicts 156.62 KN compared with a measured tension of 165.92 KN, 
while in experiment 2 the calculated tension is 162.56 KN with a measured 
165.5 KN. Although the mathematical curve has to have the same behavior 


than before (in figure 


general, to the simuls^ 

/ 


larger stretch of the v\ 
causes a bit larger p 
than in Exp. 1: 1.5 mr' 



are closer to each other, and in 
lotion for all these effects is the 
b little increase of velocity which 
larger pull through the sockets 
jf mm in exp. 2. This also means that 


a larger amount of energy was dissipated resulting in a minor measured 
tension as was observed by CP.tfatblle too, but not considered by the model, 
[31]. Figure 36, shows the wire rope pull through the sockets 13 . 


The difference in tensions also causes an over-prediction of the 
maximum acceleration along the unloading stage, as showed in figure 35. 
And due to the bigger difference in acceleration, velocities also differ being 
lesser than the experimented while loading and greater when unloading (in 
module), then the displacement, figure 33 is smaller, specially, in the 
unloading stage. Also, as with experiment 1, the calculated zero velocity 
happens almost at the same time than the experimented but does not 
coincides with the maximum tension. 


13 It is worth noting here that the elastic wave refraction in the anchorages, which prevents reflection 
to duplicate the wave amplitude, can be included in the developed equations following the 
geometrical method presented in 4.1. 
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In sum, it is clear that when a wire rope is tensioned up to its plastic 
zone and with a more changed elastic modulus and additional dissipation in 
its anchorages, then the mathematical model could predict the maximums 
but over-predict, with some percentage, tension level in both stages of the 
movement. 



sockets filled with resin used as anchorages in the experiment 2 of C. 
Stolle [31]. 

Something to note here is how the numerical simulation of C. Stolle also 
shows a greater difference with the experimental data, compared with his 
simulation of experiment 1. He also explain this with the larger permanent 
stretch of the wire rope and energy dissipation on the anchorages [31]. 
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5.2.2 . Modeling experiment 4 of C. Stolle. 

Entry variables and parameters to model the experiment 4 of C. Stolle 
are in tables 18 and 19. Calculation results are presented in figure 37, 38, 39 
and 40 and tables 20, 21,22 and 23. 


Distance Initial wire rope 

Impact velocity Mass of the car between tension (New wire rope 

anchorages without pre-tension) 

1.06300E+01 m/s 7.69300E+02 Kg 8.61100E+00 m 0.00000E+00N 

Table 18: Entry variables to model experiment 4 of C. Stolle [31]. 


Impact body 
(car) mass M c 

7.6930E+02 Kg 


Wire rope 
metallic area 

1.54508E-04 m 2 


Impact body 
velocity V 00 

1 .05700E + 01 m/r, 8.61100E+00m 


Ideal cable length 
L. 


Elastic modulus 
E 

1.04000E+11 Pa 




M= M + M hl 

_c_ cable. 

7.79875E+02 Kg 


Table 19: Parameter 
rope loading stage. 



^ nd velocity in 
steel C„ 

- 8E+03 m/s 


Wire rope mass 

M hl 

cable 

1.05746E+01 Kg 


U. 


t 


c.MAX 


84E+02 m 05 /s 7.93582E-02 Seg 
nent 4 of C. Stolle for the wire 
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7.33976E-02 Seq. 1.04000E-01 Seq. 8.62000E-02 Seq. 6.62788E-02 Seg. 

Table 20: Comparison between the maximums of the tension at the anchorages of the 
transversely impacted wire rope of the Exp. 4 of C. Stolle. 
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Table 21: Comparison of maximums displacements of the impact body, calculated with the 
measured and simulated in the Exp. 4 of C. Stolle. 

















0.05 0.1 0.15 


BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DE INFORMATICA Y COMUNICACION 



CO 


o 



5 



4-* 



c 

CD 

-S 

0 

03 

5 ns 

0 

E 

0 

CO 

O 

o 

CM 

■ — -U 

o 

o 

O 


+ 

LU 

LU 

Sis 

■o -E 

O 

LO 

LO 

o 

o 

o 

CM 

0) 0 

o 

CO 

CO 

3 

0 

CO 

CO 

cc 



CD 






> d 

~ X 
~ LU . 

0 

03 

0 

■2 <D « 

E 

CO 

> £ o 

o 

o 

o 

E=w 

+ 

LU 

1 

LU 

£ "C (j 
0 £ Z. 
QC 3 o 

O 

o 

o 

O 

03 

CO 

O) 

Ctf 0 

*cr 



T_ 

2 £ 






-4-J = 

o O 

X ■4— t 

0 

03 

0 

£ t/) 

E 

CO 

0 . 

o 

T — 

> o 

o 

o 

c > 

+ 

LU 

o 

CO 

O 

N- 

LU 

o 

o 

CO 

CO 

o 

>5 = 


T— 

!.§ 



0 



■E'TS 


03 

0 

o o 

0 

O ^ 

E 

CO 


o 

o 

CM 

o 

turn 
e Ma 

y’ M - 

+ 

LU 

*cr 

CO 

1 

LU 

CO 

1— 

0 s. 
od ~ 

CO 

CM 

o 

03 

CO 

CO 

s & 

03 

1^ 

2 0 



nS 




w 

© 


in 

d 


cn 


LO 

co 


00 


LO 

1^ 


(s/iu) w.A A);oo|3A 


O) 

co 

0) 

3 

3) 


T3 

C 

CD 

T3 

CD 

i— 

3 

CO 

CD 

CD 

E 

CD 


■a 

CD 

-t—> 

_CD 

3 

_o 

CD 

O 

>. 

■a 

o 

.a 


o 

CD 

CL 

E 

CD 


CO 

.0 

-t—> 

'o 

_o 

0 

> 

c 

i— 

3 

-t—< 
0 


_0 

"5 

-t—> 

C/) 

d 

>+- 

o 

CL 

X 


0 

E 

| LU 

II 


si 

I— 0 


Esfa ofcra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Peru. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/Hcences/by-nc-sa/2.5/pe/ 



















0.02 0.04 0.06 0.08 0.1 0.12 


BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DE INFORMATICA Y COMUNICACION 


o 


<\l 


00 VO 


<N 


CD 



c 

CD 


§ « 


_ O’ 


£ w 


.2 ^ ^ 

(C N 

CB a) c/> 

O) 

^ 1. — 

o 

O (fl h 

+ 

H -O 00 

LU 

O r- 

cd E o 

Co || 

CO 

O 

Ia» 

■'cf 

o 

3 O 
« 5 
cB 5 


CD 




T5 

CD 1 1 rH 


= O oo 

SO?? 

CD 

t) N Ol . 

o 

^ duf; 

a> uj ®E 
o 

+ 

LU 

o 

CO 

<r a> cu 

CD 

^ s. — 

CO 

o ~ o 

o 

3 »- 
.= (/) 

T“ 



>■» 

X! £ 


T5 H 


0) 00 


■*-» ■ 
cB 9? 

CD 

— ii 

T— 


O 

£ - ■ 

+ 

<55 «sr 
. D) w 

LU 

O 

CD 

s I 

"sf 

O 

< O 

T 

. +; 


X 0) 


CB . 

S o 


® tfT 
+-■ 05 


'' 

CD 

as « 

O 

c/> > ■£ 
<® = 

<D O 

O 

+ 

LU 

00 

00 

O 5 O) 

CM 

< ^ II 

"sf 

■ £ 05 

X £ 

CB CB 

o 

2 2 



C\J 

CO 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Peru. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 


7.33976E-02 Seg. 8.79130E-02 Seg. 8.66130E-02 Seg. 8.62000E-02 Seg. 

Table 23: Maximums accelerations of the impact body, calculated with the measured and simulated in the 
Exp. 2 of C. Stolle. 
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5.2.2.1 .Discussion of the modeling of experiment 4 
of C. Stolle. 

Here, only a loading stage is present, up to two strands break, this is 
why is used an elasticity modulus of 104 Gpa. Mathematical predictions of 
tension, acceleration, velocity and displacement have similar characteristics 
than those made for experiments 1 and 2: modeled acceleration is greater 
which results in a predicted velocity lesser than the measured while loading, 
although the displacement is very close up to the time when it breaks. 
Calculated tension is also greater because the model does not include the 
permanent stretch and the pull through the sockets. This is also why the 
predicted time when it breaks is less than the measured or simulated, then 
the maximum accelei 
accurately. But is the 
break what provokes tli 
rope breakage. C. Stc" 


orfc ilacement can not be predicted 
„ r he wire rope up to two strands 
model can only predicts the wire 
>ate the behavior as well, and the 
basic reason is the n< wire ro P e 9 eometr y simplified to 

an uni-dimensional string or a set of points like C. Stolle does. Figure 57 
shows how the experiment fin&Sgifr with only one strand of the 3x7 wire rope 
without breaking but stretched. 



However, a measure to take in count is the level of tension at which the 
two strands break, 171.26 KN, while the maximum tension it can withstand is 
181 KN. But a difference of ~10 KN or 5.5% is not big for an engineering 
application because the security factor used for wire ropes is always more 
than 2.5 [75], In addition, this confirm that assuming the anchorage tension 
as equal or similar to that of the impact point is very close to what reality 
shows and confirm that the wire rope inertia is not an important factor to 
consider when it is impacted with the aforementioned conditions. Although it 
does indicates a decrement in the wire rope breaking strength at the impact 
point and this is the main reason why the internal stresses model was 
developed. 
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5.3 . Validation of the model of longitudinal stresses 
and the model of the impact body’s movement 
using the simulation data of Tytko et al and P. 
Gospodarczyk. 

Both researches used a 1x7 wire rope, its characteristics are in table 24 


and figure 41: 


Wire rope type 

Number of strands 

Total number of 
wires 

Nominal diameter 

1x7 

1 

7 

15.7 mm 

Wire rope metallic 
area 

1.50000E-04 m 2 

Weight per unit 
length 

1.17200E+00 Kg/m 

Steel density 

7.80000E+03 Kg/m 3 

Elastic limit 

2.46000E+02 N 

Breaking strength 

2.79000E+05 N 

Elastic modulus of 
a wire repo&S&i'&i/ 

nqyire^tesl/u 

1.9# n 

!TV 



Table 24: Physical cha 
al and P. Gospodarczy 



3teel wire rope used by Tytko et 


5.233 mm 

Figure 41: Cross section of a 1x7 wire rope 
with 015.7 mm. 

Both works used the software LS-DYNA (also used by C. Stolle), the 


authors do not give any detail about their numerical models, but they used 
the following conditions for the wire rope: An completely straight wire rope 
(simplified), its length was exactly equal to the distance between anchorages 


which were perfectly rigid. 
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Simulation of Tytko et al is modeled first, he does not present results for 
the kinematics of the impact body. Data and parameters to take in count are 
the wire rope properties of table 24, the block impact velocity V 00 , its mass 


M c and the distance between anchorages given in tables 25 and 26 and 


showed in figure 42. 


M c =4000 Kg 
Rigid anchorages 


V 00 = 3.0 m/s 


Wire rope 1x7 


s 


±- 





Figure 42: ldealizat : o^prT 
length equal to the ,r- 
Voo, is the velocity 
origin of coordinates 
Tytko et al, pres< 

rope breaks under the 


"^^f Tytko et al. L is the wire rope 
**C, is the impact body mass, 
‘l§prdinated axis and O, is the 
wire rope. 

3 showing how his modeled wire 
cted. 


\ integration [ 
min=-0.678102, at elem# 2952 
max=1.95835, at elem# 2970 



Figure 43: Transverse impact over a wire rope and its breakage 
simulated by Tytko et al using LS-DYNA. 


Below, the developed mathematical model is applied for the wire rope 
tension at the anchorages (or longitudinal) and compared with the simulation 
results obtained by Tytko et al. The work of Gospodarczyk requires another 
analysis this is why is given in the next section. 
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5.3.1 . Modeling the simulation of Tytko et al. 

Entry variables and parameters are in tables 25 and 26. Calculation 
results are presented in figure 44 and tables 27 [29], 


Impact velocity 

Block mass 

Distance between 
anchorages 

Wire rope initial 
tension 

3.00000E+00 m/s 

4.00000E+03 Kg 

I.OOOOOE+OOm 

0.00000E+00 N 


Table 25: Entry variables to model the simulation of Tytko et al [29], 


Impact body 
velocity V 00 


Impact body 
(block) mass M c 

4.00000E+03 Kg 

Wire rope 
metallic area 

1.50000E-04 m 2 
M= M + M hl 
4.00117E+03 Kg 
Table 26: Parameter 
rope loading stage up 


Ideal cable length Elastic modulus 
L. E 



1.00000E+00 m 


1.95000E+11 Pa 


^■■nd velocity in 
steel C M 

POE+03 m/s 

U„ 


Wire rope mass 

M M 

cable 

1.17000E+00 Kg 


01E+02 m 0 5 /s 


tion of Tytko et al for the wire 
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1.99049E-02Seg. 1.98825E-02Seg. 2.79000E+02 N 2.89630E+05N 8.74344E-01 mm/mm/s 

2.46000E+05 N 2.46000E+05 N 2.12365E-02 Seg. 3.80749E-02 Seg. 2.12153E-02 Seg. 

Table 27: Times to reach the elastic limit with the calculated and simulated breaking strength for the simulation of Tyko et al. 
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5.3.1.1 .Discussion of the modeling of Tytko et al’s 
simulation. 

In this case only calculation of tension is necessary and then its 
comparison with the simulated data of Tytko et al. The closer approximation 
is due to the more ideal conditions imposed in Tytko et al’s numerical model: 
the anchorages are truly rigid and the wire rope length coincides with the 
distance between them, also the elastic modulus is 195 Gpa, which is very 
close to that of a steel bar: 207 Gpa, this is, is close to a constant. Hence and 
for these reasons, the mathematical model predicts accurately the simulation 
until the elastic limit of reference, 246 KN, where the slope of the numerical 


curve begins to decrease due to the plastic behavior of the wire rope 


until its 


breaks at 289.63 KN. 
which does not con 


other hand, in the 
developments, the ela 
the simulation. 



Ily greater than its static strength 
high deformation ratios. In the 
an be included, in posterior 
Stic zone and approximate better 


m 

An important observation here i§ T#w, even if the impact body is not a 
material point or does not have a narrow cylindrical front, instead it has a 
plane front, the mathematics is able to approximate the numerical result 
mainly because the large differences between the mass of the block and the 
wire rope. 
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5.3.2 . Modeling the simulation of Gospodarczyk. 

Entry variables to model this case are in table 28. 


Impact velocity 

Parallelepiped 

mass 

Distance between 
anchorages 

Wire rope initial 
tension 

7.50000E+01 m/s 

1.80000E+04 Kg 

7.50000E+00 m 

0.00000E+00 N 

Wire rope additional 
length 

Mass-wire rope 
contact shape 

Initial wire rope 
shape 

Number of wire 
ropes in the base 


Plane front of a 



O.OOOOOE+OOm 

parallelepiped (See 
figure). 

Straight 

420 


Table 28: Entry variables to model the simulation of P. Gospodarczyk [33], 



Figure 45 shows the initial position, an instant before the impact 
between the 18000 Kg p' J c 1.5m at 75 m/s and an artificial 

base (barrier) that ^ % ( w 'th a non specified small 

thickness, the author J muc h) with 420 1x7 wire ropes 

in a criss-cross patte s ’ x ' eve ' s 6*(2*13+4*11). The 

parallelepiped contact SvDli 

1 


6 levels of* 
wire ropes 


07.5m 


Figure 45: Initial position of the parallelepiped and base made of 
wire ropes an instant before the impact simulated by P. 
Gospodarczyk. The base has diameter of 07.5m and is used in 
deep shaft mines. 1 is the parallelepided, 2 is the set of wire 
ropes and 3 is a sheet metal [33], 


The idealization begins supposing that only the wire ropes directly 


under the parallelepiped suffer the impact and that the sheet metal strength 
is negligible. This idea is shown in figure 46: 
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13 1x7 wire ropes 



(13 + 8) x 6 = 126 1x7 
wire ropes under the 
parallelepiped. 


13 x 6 1x7 wire ropes 
Under the side of 2.5m. 


Figure 46: First simulation of an impacted 


artificial base u 
Gospodarczyk. 

Then, to be abU 

done, this way: the m 


mines and simulated by 
k % 

a second idealization has to be 
ded equally between all the wire 
. 86 Kg per wire rope. And now 


ropes directly under it 
this is the mass of the material point tq&e used into the models impacting a 
1x7 wire rope at 75 m/s as illustrated in figure 47. 


V Q0 = 75.0 m/s 


Rigid anchorages 


"M c =18000 Kg /126 wire ropes 
M c =142.86 Kg 

Wire rope 1x7. 

—► 
X 


ly 

L = 7.5 m 


Figure 47: Second idealization of the simulation of an impacted artificial 
base used in deep shaft mines and simulated by Gospodarczyk. L, is the 
distance between anchorages that coincides with the diameter of the 
base, M c , is the mass of the block, Voo, is the impact velocity of the block, 
x is the coordinated axis and O, is the origin of coordinates and the middle 
point of the wire rope. 
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Then, as in the rest of cases, equation (39) is used, to calculate the 
tension. But here equation (29) of the model of small displacements is used 
to calculate the velocity of the block during the impact. All the reasons about 
why the equation of large displacements is not used are given in the 
discussion of this case in 5.3.2.1. Entry variables and parameters are in 
tables 28 and 29, calculation results are in figures 48, 49 and table 30. 


Impact body 
mass M c 

Velocity of impact 

v 00 

Wire rope 
idealized length L. 

Elastic modulus 

E 

1.42860E+02 Kg 

7.50000E+01 m/s 

7.50000E+00 m 

1.95000E+11 Pa 

Wire rope 
metallic area 

Wire rope steel 
density 

Sound velocity C lf 

Wire rope mass 

M hl 

cable 

1.50000E-04 m 2 

7.80000E+03 Kg/m 3 

5.00000E+03 m/s 

8.77500E+00 Kg 

Mass in 

movement M c + 

^cable. 


u 

% 2 

t c .MAx@278.6 KN 

1.51635E+02 Kg 

8.5 MSjmrkm 

02E+02 m 0 5 /s 

7.33144E-03 Seg 


tc.MAx@278.75 KN 

7.33462E-03 Seg 


Table 29: Parameter: 
ropes in an artificial b 


to moc 
of a dee 


on of Gosodarczyk for the wire 
ne during the loading stage. 
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Table 30: Final velocity calculated and simulated of the 18 TN block after 
passing through the artificial base simulated by Gospodarczyk. 
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5.3.2.1 .Discussion of the modeling of 
Gospodarczyk’s simulation. 

As here is used the same 1x7 wire rope as in the case of the work of 
Tytko et al, then the same elastic modulus is chosen, 195 Gpa, and also the 
same breaking strength (static) of 279 KN. The sole objective in this case is 
to find the final velocity of the parallelepiped after going through the entiry 
base breaking the sheet metal and the wire ropes. First, the time when the 
wire ropes break is calculated and with this is calculated the block velocity 
reached just after it breaks the last wire rope of the base. But to calculate this 
velocity is not a simply use of equations. Here the impact velocity of 75 m/s 


have to 


and the parallelepiped dimensions compared to the wire ropes length 
be considered before eloc% corresponds in the third scale of 


the Weihrauch chara 
present and up to th 
simulation and with a' 

\ \ v 

covers a considerable 


'.Stic deformation is going to be 
as is observed in Tytko et al’s 
!>sion. Also, the block plane front 
s length of 7.5m: the 1.5m side 


covers a 20% and the 2 5m side a 55.5%. Therefore during the impact the 
wire rope points directly unJi' the block only have deformations in an 
horizontal direction and its slope is going to be always zero or with a direction 
parallel to the plane surface of the rigid parallelepiped. The difference is how 
an ideal elastic string deforms when is impacted by a block with a shape or 
by a material point, in this situation the local inertia of the wire rope and its 
change driven by the stress waves provokes that the slopes of every point 
around the impact zone to be greater than the ones around the anchorages 
as is observed in figures 11 and 14, also N. Cristescu shows this on his 
numerical solution of the impact over a string with one end free and C. Stolle 
on the photographs of his experiments [31], [37], [53], But these differences 
are not present when the parallelepiped impact the base, there, the points in 
the impacted zone stretches but do not contribute to the vertical 
displacement of the wire rope. 
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In sum, this is the answer to why is used the model of small 
displacements (light blue delimiter string) to approximate the results of the 
simulation with a predicted final velocity of 51.29 m/s very close to a 
simulated result of 48 m/s. While if the large displacements equations (62) 
and (65) are used then the final velocity predicted is slightly different to the 
impact velocity. The error in the approximation with the small displacements 
model happens because it does not include the sheet metal, any plastic 
deformation or any steel hardening by the high deformation ratios. Figure 50 
shows the simulation result. 



Figure 50: Final result of the simulation made by Gospodarczyk 
with the block decelerated from 75 m/s to 48 m/s [33]. 
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5.4. Qualitative comparison between the internal 
stresses model application and experimental 
and simulated observations. 

In this part, the experimental observations are the ones of C. Stolle and 
from the simulation only of P. Gospodarczyk’s work [31], [33], Simulation of 
Tytko et al is not considered because they do not give the impact block 
measures. 


Here the model of internal stresses by the contact during the 
transverse impact over a wire rope is used. All the variable and parameters 
are used from the experiments and simulations specified above. Tension in 


the impact zone is assume^^Qlii^ reasons presented before, equal to the 

Sr*' . 

in this case only the maximum 

[V \ 

^contact force F Mdx predicted or 


tension in the anchor 


value is used T 


Max 


measured (maximum 
the model, a circular i 
equivalent to its meta 



by mass). Then, in order to use 
le cross section of the wire rope 
le cross sectional areas of each 


of its constituent wires, /^maximum calculated acceleration is found from the 
maximum prediction during the loading or unloading using equations (69) and 
(70), the equivalent diameter is found with equation (72), the internal stresses 
along the transverse z axis of the equivalent perpendicular cylinders is 
found with equations (76), (77) and (78). Internal stresses for a plane rigid 
body impacting with a wire rope are found with equations (80), (81), (82) and 
(83). Finally an equivalent stress which combines the internal stresses in the 
contact zone with the longitudinal stresses are calculated using the Von 
Mises criterion with equation (84). 
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5.4.1 . Application of the internal stresses model to 
the experiment 1 of C. Stolle. 

Conditions and parameters of C. Stolle experiments are in 5.1. The 
experimental car he used had a front with a cylinder of 076mm made of A36 
steel as figure 51 shows [31]. 



Figure 51: Cylindrical impact front used in the experimental car 
of C. Stolle [31], 


Then the maximum contact force and tension for each case have to be 
calculated using the indicated equations. The 3x7 wire rope characteristics 
are in table 3 and figure 25. The same elasticity modulus used before to 
model the loading stage is used. A breaking strength o^ MAX is calculated 
with the breaking tension indicated by C. Stolle, 181 KN divided by the wire 
rope metallic area. Entry variables from experiment 1 of C. Stolle are in table 
31, calculation results are in figure 52 and table 32. 
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Wire rope metallic area 

1.5451E-04 m 2 

Wire rope equivalent diameter 

1.4026E-02 m 

Cylindrical impact body diameter 

7.6000E-02 m 2 

Wire rope Poisson ratio 

3.0000E-01 

Impact body Poisson ratio 

3.0000E-01 

Wire rope elastic modulus 

1.0400E+11 Pa 

A36 steel elastic modulus 

2.0700E+11 Pa 

Maximum tension measured: T c , ... 

Expl. Max. 

6.7395E+09 N 

Maximum tension calculated: T . , ... 

mod. Mat. Expl. Max 

6.6538E+09 N 


Table 31: Entry variables to apply the model of internal stresses 
to the experiment 1 of C. Stolle. 

Wire rope tension: Measured (Med.) and 
calculated (Cal.) 


T m . =165.2 

T :a,c. =156.62 KN 

Equivalent diameter of the cylinde 


ntact 

1.6325E-02 

1.6325E-02 m 

Radio of the contact area a^f^ r ' 



2.2719E-03 

2.2430E-03 m 

Contact pressure P 0 . ^ 



6.7395E+09 

6.6538E+09 

Maximum contact force Mol | 

F . c =M*a„„. _ , J/Wi. 

max.Exp2 Max.Exp2 jXjm 


ik 

1m 

7.2855E+04 

7.0113E+04 N 

Maximum longitudial stre/ 

n jMr/1 

1.0692E+09 

1.0137E+09 Pa 

Superficial Von- Mises sti 


pyo) 

2.0979E+09 

2.0365E+09 Pa 

Maximum internal Von Mi 

nu 

4.7992E+09 

4.7082E+09 Pa 

Minimum internal Von Mise«;stress 

•fr- V.Min 


1.0869E+09 

1.0307E+09 Pa 

Breaking strength ct rm4x 



1.1715E+09 

1.1715E+09 Pa 

% of the wire rope with stresses greater than the 
breaking strength over the axis Z of its cross 
sectional area. 

43.00 % 

35.00 % 


Table 32: Calculation results of the internal stresses model using data of the 
experiment 1 of C. Stolle. 
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5.4.2 . Application of the internal stresses model to 
the experiment 2 of C. Stolle. 

Entry variables from the experiment 2 of C. Stolle are in table 33, 
calculation results are in figure 53 and table 34. 


Wire rope metallic area 

1.5451E-04 m 2 

Wire rope equivalent diameter 

1.4026E-02 m 

Cylindrical impact body diameter 

7.6000E-02 m 2 

Wire rope Poisson ratio 

3.0000E-01 

Cylindrical impact body Poisson ratio 

3.0000E-01 

Wire rope elastic modulus 

1.0400E+11 Pa 

A36 steel elastic modulus 

2.0700E+11 Pa 

Measured maximum tension: T c , u . 

Exp2. Max. 

1.6550E+05 N 

Calculated maximum tension: T 

^ i\ir|r> Mat Ex P 2 - Max 

1.6256E+05 N 


Table 33: Entry 
to the experimen 

/f 


model of internal stresses 


Wire rope tension: Measi MMfak 
calculated (Cal.) 

lli V 


=165.5 =162.56 KN 

Equivalent diameter of tf 



1.6325E-02 1.6325E-02 m 

Radio of the contact aree 

2.2885E-03 2.2849E-03m 


Maximum measured contaci force in Exp.2 
F =M*a 

max.Exp2 Max.Exp2 


7.4467E+04 7.4113E+04N 


Maximum longitudinal stress o L , 


1.0711E+09 1.0521E+09 Pa 


Superficial Von Mises stress a, 


S.Mtn 


2.1084E+09 2.0907E+09 Pa 


Maximum Von- Mises stress ct, 


V.Max 


4.8347E+09 4.8157E+09 Pa 


Maximum Von- Mises stress cr, 


V.Min 


1.0893E+09 1.0702E+09 Pa 


Wire rope breaking strength a 


R.Max 


1.1715E+09 1.1715E+09 Pa 


% of the wire rope with stresses greater than the 
breaking strength over the axis Z of its cross 
sectional area. 


44.00 % 


40.60 % 


Table 34: Calculation results of the internal stresses model using data of the 
experiment 2 of C. Stolle. 
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Figure 53: Results of the model of the internal stresses (maximums) around the impacted zone of the wire rope used in the 
experiment 2 of C. Stolle. With two curves one using the maximum calculated tension Tmod. m at.Expi. Max , and the other with the 

measured tension TExpi.Mdx and compared with the wire rope breaking strength a R . M a X . 
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5.4.2.1 .Discussion of the application of the internal 
stresses model to the experiments 1 and 2 of 
C. Stolle. 


Both predictions are similar. Using the measured data, results says that 
a 43% in exp. 1 and a 44% in Exp. 2 of the wire rope cross section in the 
impacted zone would to be broken. Whereas by using the calculated data the 


percentages decrease to 35% and 40% respectively. Although, nothing of 
these predictions happened, as are observed in figures 54 and 55. The wire 
ropes ended curved around the impact zone and none strand (which are a 
33.33% of the wire rope each one) broke. Thus, the model is not able to 



predict the experimental results in these cases. 


Figure 55: Permanent 
deformation of the wire 
rope around the impact 
zone after the Exp. 1 of 
C. Stolle [31]. 


Figure 54: Permanent 
deformation of the wire 
rope around the impact 
zone after the Exp. 2 of 
C. Stolle [31]. 


This error happens because the model does not take in count the 
rearrangement of wires and strands when are crushed by a transverse force. 
Stresses does not distribute as in a continuous cross section, its shape 
changes continuously while being impacted. The 3x7 wire rope has an 
arrangement showed in figure 25, each one of the three strands is held by 
the others only in two lines along a spiral, with no restriction to its lateral 
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movement 14 . This implies that there is not a continuous path where stresses 
can be transmitted all over the cross section as could happen to a larger 
extent in a 18x7 wire rope, for example, or in any other wire rope with more 
strands without a synthetic or natural core. 


5.4.3 . Application of the internal stresses model to 
the experiment 4 of C. Stolle. 

Entry variables to model the internal stresses are in table 35, calculation 
results are presented in figure 56 and table 36. As in experiment 4 of C. 


Stolle the wire rope breaks then only the loading elastic modulus is used and 


the same breaking strength are used as before in the other cases. In table 
37, are the calculations usina two aiiernative tensions: one equal to the 


maximum measured Jp 
maximum accelerate j 
strength changes ther„ 


v % higher. And are used with two 
dingly (because the breaking 
o changes). 


Wire rope equiva 

?n / / 

1.4026E-02 

m 2 

Cylindrical impac. i^ody d. 


7.6000E-02 

m 2 

Wire rope Poisson iiftffl 


3.0000E-01 


Cylindrical impact bodyPSJiiasbl 

iT*fatio 

3.0000E-01 


Wire rope elastic modulus 


1.0400E+11 

Pa 

A36 steel elastic modulus 


2.0700E+11 

Pa 

Maximum measured tension in Exp4: T c . 

r Exp4. Max. 

1.7126E+05 

N 

Maximum breaking tension measured statically 
by C. Stolle: T . _ . ... 

J mod. Mat.Exp4. Max 

1.8100E+05 

N 

Maximum alternative tension: T 

max. Alt-1 

1.7490E+05 

N 

Maximum calculated tension: T 

max. Alt-2 

1.7126E+05 

N 


Table 35: Entry variables to apply the model of internal stresses 
to the experiment 4 of C. Stolle. 


14 A transverse force over a wire rope can also happens in other cases, for example a flattening can be 
caused by bending the rope severely over the rim of a sheave or any sharp object [76]. 
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Maximum wire rope tensions: Measured Exp.4 
(Med.) and calculated (max. Strength) (Cal.) 

T M ed =171-26 

T ca,c. = 181 KN 

Equivalent diameter of the cylinders in contact 

1.6325E-02 

1.6325E-02 m 

Radio of the contact area: a. 

2.3542E-03 

2.0545E-03 m 

Contact pressure P 0 . 

6.9837E+09 

6.0947E+09 Pa 

Maximum measured contact force Exp.4 

F . _ ,=M*a M . _ , 

max.Exp2 Max.Exp2 

1.7126E+05 

1.8100E+05N 

Maximum longitudinal stress o L M4x 

1.1084E+09 

1.1715E+09 Pa 

Superficial Von- Mises stress tJ SMin 

2.1743E+09 

2.0703E+09 Pa 

Maximum Von- Mises stress 

4.9478E+09 

4.4811E+09 Pa 

Minimum Von- Mises stress 

V.Min 

1.1282E+09 

1.1845E+09 Pa 

Wire rope breaking strength ct r MSx 

1.1715E+09 

1.1715E+09 Pa 

% of the wire rope with stresses greater than the 
breaking strength over the axis Z of its cross 
sectional area. 

58.00 % 

100.00 % 


Table 36: Calculation r 
experiment 4 of C. Stol 

/#' 


Is of theyinterf 


itresses model using data of the 


Wire rope maximum 

T =174.9 T 

max. Alt-1. r 

rax. Alt-2. 171.26 KN 

Equivalent diameter of t 

1.6325E-02 

1.6325E-02 m 

Radio fo the contact ar 

2.0199E-03 

1.9988E-03 m 

Contact pressure P 0 . 

5.9921E+09 

5.9293E+09 Pa 

Maximum calculated contaC^g^jE^^S^ 

F =M*a 

max.CalcExp2 Max.Calc.Exp2 

5.1205E+04 

4.9612E+04 N 

Maximum longitudinal stress o L M . x 

1.1320E+09 

1.1084E+09 Pa 

Superficial Von- Mises stress cr SM . n 

2.0185E+09 

1.9873E+09 Pa 

Maximum Von- Mises stress a..... 

V.Max 

4.3916E+09 

4.3369E+09 Pa 

Miminum Von- Mises stress 

V.Min 

1.1444E+09 

1.1204E+09 Pa 

Wire rope breaking strength a^ Max 

1.1715E+09 

1.1715E+09 Pa 

% of the wire rope with stresses greater than 
the breaking strength over the axis Z of its 
cross sectional area. 

58.00 % 

45.00 % 


Table 37: Calculation results of the internal stresses model using data of the 
experiment 4 of C. Stolle using two alternative maximum tensions. 
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Deep from the origin of the Z axis (% of the wire rope diameter) 

Figure 56: Results of the model of the internal stresses (maximums) around the impacted zone of the wire rope used in the 
experiment 4 of C. Stolle. With two curves one using the maximum calculated tension Tmod. m at. Exp i. Max , and the other with the 

measured tension Tex P i. Max and compared with the wire rope breaking strength o R . M a X . 


















BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DE INFORMATICA Y COMUNICACION 


5.4.3.1 .Discussion of the application of the internal 
stresses model to the experiment 4 of C. 
Stoile. 

Figure 56 shows how, when the wire rope tension at its anchorages is 
171.26 KN (maximum measured tension) then the model predicts that a 58% 
of the cross section has to be broken which almost coincide with the 66.66% 
that correspond to two strands broken as happened in the experiment 4 of C. 
Stoile. In the other hand if the maximum breaking strength is used, 181 KN 
all the wire rope have to be broken an this coincides with the supposition of 
the analytical model. And if the alternative maximum tensions are used with 
their correspondent calculated maximum accelerations then by using a 



|sured) then only a 45% would be 


>d then a 58% broken section is 
\ :ed using the experimental data. 


^^scoincide exactly with the reality 
@ ken section, by using a tension of 
iv referential and predicted breaking 


strength of the impacted wire 


However, as calculations of experiments 1 and 2 are far from reality 
then the internal stresses model is inconsistent. And an explanation to why 
there is a coincidence with the results of experiment 4 is: first the 
rearrangement of wires and strands is done during the first milliseconds, 
next, a reconfigured cross section bears the impact force with small changes 
but as the impact velocity is 10.57 m/s, which is up to 56.6 % greater than the 
impact velocities of experiments 1 and 2, and with the tension reaching a 
maximum of 165.5 KN, then the breakage is inevitable. The particularity here 
is that the breaking tension was 171.26 KN which is not the maximum static 
measured of 181 KN and is also slightly bigger than the maximum tension in 
the other mentioned experiments. But at the end is enough to coincide with 
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the prediction of the internal stresses model. Only because the breaking 
strength is close to the achieved with a much less velocity. 

Figure 57 shows the final condition of the wire rope after finishing the 
experiment 4 of C. Stolle. Only one strand of the 3x7 wire rope stand 
unbroken and locally bent. 
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5.4.4 . Application of the internal stresses model to 
the Gospodarczyk’s simulation. 

Characteristics of the wire ropes and the simulation conditions are in 
5.3.2. The impact solid has a plane front of 1.5m by 2.5m and falls over wire 
ropes of 7.5m long. Here, as in the simulation, the solid is considered rigid 
and its mass is divide between the 126 wire ropes (see 5.3.2.1) [33], To use 
the model is only necessary the final contact force and the reached tension. 
Entry variables are given in table 38, calculation results are presented in 
figure 58 and table 39. 


Wire rope metallic area 

1.5000E-04 

m 2 

Wire rope equivalent diameter 

1.3820E-02 

m 

First side of the rigid parriift§l&#]^ed n-1. 

2.5000E+00 

m 

Second side of th ,^cl -2. 

1.5000E+00 

m 

Wire rope Poisson *7 ;y 

3.0000E-01 


Wire rope elastic ' 

1.9500E+11 

Pa 

Maximum calcula Matcosp" 1 

2.7875E+05 

N 

Maximum calcula ....__-2 

2.7856E+05 

N 


Table 38: Data km Gwupwuui Ks£—y ^ simulation needed to use 
the model of internal 'tresses. 


Calculated breaking tension for: side of 1.5m and 2.5m. 

Tca.c.sm =278.56 

T calc. 2.5m =278.75 KN 

Half width of the contact area b. 

2.1989E-04 

1.7033E-04 m 

Contact pressure P 0 

5.4554E+08 

1.1739E+09 Pa 

Calculated maximum contact force 

F =M*a 

max.Gosp. Max.Gosp. 

4.7108E+05 

4.7112E+05 N 

Maximum longitudinal stress ct l m4x 

1.8583E+09 

1.8571E+09 Pa 

Superficial Von- Mises stress a s Min 

2.0765E+09 

2.3266E+09 Pa 

Maximum de Von- Mises stress ct„ ... 

2.0765E+09 

2.3266E+09 Pa 

Minimum de Von- Mises stress o\,... 

V.Min 

1.8600E+09 

1.8600E+09 Pa 

Maximum breaking strength ct r MAx 

1.8600E+09 

1.8600E+09 Pa 

% of the wire rope with stresses greater than the 
breaking strength over the axis Z of its cross 
sectional area. 

100.00 % 

100.00 % 


Table 39: Calculation results of the internal stresses model using data of 
Gospodarczyk’s simulation. 
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Figure 58: Results of the model of the internal stresses (maximums) in the impacted zone of the wire ropes of the artificial base 
for deep shaft mines simulated by Gospodarczyk using two maximum calculated tensions Tmod. m at.Expi. Max , one per each side of 

the parallelepiped that impact and compared with the wire rope breaking strength o R . M a X .. 
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5.4.4.1 .Discussion of the application of the internal 
stresses model to the simulation of 
Gospodarczyk. 

Prediction using the models indicates that the wire rope capacity due to 
a concentration of stresses at the impact (or contact) zone is not 
considerable. Statically measured the standard breaking strength of the wire 
rope is 279 KN then for the wire ropes under the 1.5m side of the 
parallelepiped the calculated breaking stress is 278.75 KN and for the ones 
under the 2.5m side is 278.56 KN. Why this happens has its base in the wire 
construction and properties. The arrangement of six wires around a central 
one prevents, to a larger extent compared with type 3x7 (for example), its 



ach other in contact. This kind of 

m. 

cdulus of 195 Gpa which is very 


3 it does not change very much 


and new wire rope. And because 


rp# atively, what happens around the 


impact zone of the im dj* the artificial base simulated by 


Gospodarczyk. 
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6 . Conclusions and recommendations. 

Transverse impact over steel wire ropes was investigated by building 
new mathematical models using mechanic’s general laws and then validating 
them with previous experimental and numerical works. 

The general objective of determining stresses along a wire rope at its 
middle point was partially achieved: 


Two main models has been established: one of the reaction force at the 
ends or anchorages of the wire rope and another of the kinematics of 
the impact body. Both can be used either the wire rope have small or 
large displacements. A third model was done to analyze the decrement 
in strength due to^Jp'contaat strides that arise in the impact zone of 
the wire rope bug? / jjjff I k 'diet observed reality, although, is 


possible to imf 
properties The < 



. more appropriately wire rope 
n predict with more accuracy if 
the anchorages %$$$ /ire rope elasticity modules are 

constants or if it is c iaoiiizea ana is made only of steel and round wires. 

After assuming the wire rope as a cylindrical bar in contact with a 
deformable cylinder and a rigid parallelepiped, the internal stresses 
model based on a Hertzian contact stress theory was used to find the 
stresses around the contact zone of the impacted wire rope. But after 
comparing its results with experimental and numerical works, the model 
showed to be neither inconclusive nor appropriate basically because 
the parameters used like the Poisson ratio, do not represent or take in 
count all the behavior of the wire rope parts as its relative 
displacements or friction while being under the impact force. 


✓ Results of modeling the behavior of a simulated base made of wire 
ropes impacted on its center by a parallelepiped indicate that the 
models can serve as a basis to develop mathematical models of wire 
rope nets like the ones used in rock fall or landslide barriers. But more 


111 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Peru. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 






BIBLIOTECA DIGITAL - DIRECCION DE SISTEMAS DE INFORMATICA Y COMUNICACION 


experiments are needed in order to validate and know in deep the 
involved physics. 


6.1 . Conclusions of the validation of the model of 
longitudinal stresses and the movement of the 
impact body. 


Taking in count the validation using the experiments of C. Stolle and the 
simulation of Tytko et al, then the mathematical models can predict the 
anchorages tension of the wire rope and the movement of the body that 
impacts transversely it on its middle point. This happens as long as the 
anchorages are almost rigid and the wire rope modulus of elasticity is kept 
close to a constant during ;he ioadiilg^and unloading stages or if the wire 
rope is stabilized befo ^ 3e the wire used b Y C. Stolle and 

Tytko et al are of diffe respectively, then, at least for 

wire ropes with no na .'yAMi I*'' hj , helically stranded round profile 

wires and with the fore jstablished models are valid. 


And also because^experimeni 4 ot Stolle and the simulation of Tytko 
et al showed that the wire rope Liking point under impact is around the 
impact zone and with at a tension close to its static breaking strength, this 
confirm the assumption that the anchorages tension is similar to the level of 
tension along the wire rope (under the established conditions). 


And regarding the modeling of the simulation of Gospodarczyk, it is not 
possible to say that the models can model the impact over bases made of 
wire ropes. Although, it is clear that the models approximate the block final 
velocity because the particular conditions of the simulation and specifically 
the near constant properties of the simulated 1x7 wire rope. 
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6.2. Conclusions about the application of the 
internal stresses model. 


Since predictions of observations in experiments 1 and 2 of C. Stolle did 
not coincide with the predictions done with this model, then it is not suitable 
for these cases. But this does not mean that is not appropriate for others 
because, basically, influence of the selected parameters was not investigated 
and only were chosen following the referenced work. An indication of this is 
the prediction done using data from experiment 4 of C. Stolle and the 
simulation of Gospodarczyk, this showed how the calculations coincided if 
the wire rope cross section does not change significantly during the impact. 



Despite the above, it is 
should have a Poisso 
of elasticity to be i 
can be better su 
Hertzian contact 
than 3 strands and 


II metallic steel wire rope, in general, 
3 and also a transverse modulus 
, an orthotropic approximation 
say that a model based in a 
and wire ropes with more 
parts. 


Hence, more expefo$£gntal test a;e needed in order to develop or 


improve the developed model and be able to predict the additional stresses 
that arise in the wire rope due to the contact with the impact body. 


6.3 . Summary of results. 

The next has been done: 

✓ By idealizing the wire rope as a string and modeling its transverse 
impact it was possible to calculate the tension or reaction force at the 
ends of the string or anchorages of the wire rope, that resulted to be 
similar to the tension at the impact point (middle point). 

✓ Also, by idealizing the impact body as a material point and simplifying 
its interaction with the wire rope, it was possible to approximate its 
kinematics. 
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✓ Validation of the mathematical models was carried out using 
parameters and results of previous numerical and simulation works that 
did not considered analytical approaches. This way, was showed that 
the model of the longitudinal stress and of the kinematics of the impact 
body can be used to predict the behavior of a transversely impacted 
wire rope. 


✓ As part of the mathematical models construction a so called method of 
“delimiter 15 movements” was devised in order to have a path to follow 
when analyzing the movement of a transversely impacted string and 
then translated into a geometrical reference that could be translated, 
again, into mathematiG^P^Wfej/Qiathematical models of tension at the 


anchorages of t 
internal stressed 


✓ The aforementk 
conjunction with 



OUI I 


'oJcs of the impact body or of the 
he of the bodies. 

portance as it can be used in 
aracteristics, with this can be 
Mi2nsJor a string transversely impacted 


approximated the<&ord< 

7rv , 

by a material point &j©<?^plsf!e and then the strain, tension and 
kinematics of every point of the wire rope can be found by solving the 
hyperbolic partial differential equations. 


✓ A similar method was developed to model the unloading of an ideal 
impacted string, with which was demonstrated that this process is 
symmetric with respect to time of a loading process and expressed as 
the symmetric function (of time) of this process respectively. 


In sum, an engineering tool has been obtained with which is possible to 
design protection barriers and provide better and more solutions to vial 
accidents, fall of objects in mines, landslides and rockfalls. 


15 Bordering could be a better translation to “delimitante” 
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6.4 . Suggestions for future works. 

✓ In any other development or improvement to be done to the developed 
models, it should have its base in the philosophy of change, otherwise 
it could be another algorithm of solution and would not provide an 
insight about the phenomenon of transverse impact over wire ropes. 


✓ 


✓ 


✓ 


In the same way that here was developed and used the method of 
“delimiter movements” for an ideal string transversely impacted at its 
middle point, it is possible to apply the same method to model oblique 
impact and points others than the center. But before the cases with 
constant velocity must be studied, specially from this reference [35], 



When designing 
defined by the e 
how high are the 


A way to consic 
the end of the 


% iU ie security factor is going to be 
to dissipate the impact, and 
contact in the wire rope. 


es like the energy absorbers at 
a new geometrical analysis as 


presented in figure i7 and include She effect of refraction or a reflection 
that does not duplicate tension in equation (22). 


✓ Another suggestion, if the developed models are going to be used for 
design, then the wire rope answer have to be kept elastic. It is true that 
a software like LS-DYNA or another can be used but there is the work 
to find out what numerical algorithm and parameters are appropriate to 
be used in this kind of phenomena as the published works do not 
provide any code or the totality of the parameters they used. And if 
experiments are going to be done then the standard ASTM A741 -11 
(2016) have to be considered, there are suggested 3 types of wire 
ropes, but are generic denominations of various available constructions 
and this makes difficult to use only experimentation to design structures 
under impact like protection barriers. 
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8 . Annexes. 

8.1 . Annex 1. 

The aim here is to model the impact between a material point and a 
string from the first contact instant and before any wave reflects at its rigid 
ends. For this purpose figure 59 is presented: 



00 



Aproximated model 
with triangular shape 


\A/iro mno irloalivation 



impacted wire rope with and 


Ration with a triangular shape 
) es reflects at its rigid ends. M, 
actual velocity; O, origin of axis 
ige; F, represents the reaction 


force; Cm, is the longitudinal stress wave speed, X*, is the transverse 
wave speed and a is the angle ^iSfi^een axis x and the triangular profile 
sides. 

Here is represented an ideal string during its first instants of movement, 
just after it is impacted by a material point with mass M and initial velocity 
V 00 . The transverse wave has not reached the ends yet and neither the 
longitudinal wave. The curved line is the idealization while the triangular 
profile with a angles is its approximation. With this is possible to formulate 
mathematical relations as follows: 

First, the reaction force is related with the string tension with: 


F=2 Tsin a 


(85) 


But tension T can be expressed with the deformation of the triangular 
string: 
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X 

And then related with acceleration: 


( 86 ) 


F=M 


tiyo 

d.t 2 


(87) 


The extreme point x* of the triangle can be found with the transverse 
wave velocity: 


x * — 


f 


T ,* 


A cable P 


( 88 ) 


In equation (88), B , is^ge^^nt and T is the string tension and in a 


non linear relation: 


And replacing (8 



(89) 


Jtion for x* is obtained: 


**~2*%lS«^^*V 2 (90) 

Then, using equations (86) and (87) results, again, in an Endem-Fowler 

equation which models the movement of the material point: 


^ -»5/4 A 7 / 4 ~ p3/4 1 

a 1 y 0_2 A pE 


y _\ 1I2 =b U U ' 2 

r B | 13 r B| 


(91) 


d t t 2 M 

Now, by considering Yu. A. Dem’yanov’s supposition as in 4.1, the 
material point displacement has to be (like equation (11)) [66]: 

y{S 0 =0,t) = y 0 (t) = B 1 t n ;B l >0 (92) 

If this relation is replaced in equation (91), a solution is found when 


«= 2 Y B,=[BJ 2) 2 : 


y(S 0 =0,t) = (B 13 /2) 2 r 2 


(93) 
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Up to here everything seems to follow a logic path even all the 
approximations and simplifications done, but in coefficient B l3 of equation 
(91) the material point mass is dividing, indicating that while it increases then 
the material point will have less and less acceleration and as a consequence 
the reaction force F will decrease. However this has not any physical sense 
because if the body, keeping constant the velocity of impact, increase its 
mass then its momentum will increase too and after it contact an string with 
much more less mass and without bending strength, its points will accelerate 
more and with a bigger reaction force. In few words if the body increase its 
mass, it will sweep faster the string when impacting. 


An alternative is to u^srift^Btly the general law of movement of Euler or 


general law of line 

/ 4 * 

bidimensional space, t 



m unidimensional body in a 
ie movement equation is: 


= MV 


(94) 


This equation says: the impact body has a new velocity lesser than its 
impact velocity and with a new momentum that added to the string 
momentum is equal to the initial momentum of the impact body. However, 
this has to represent a wire rope being impacted and its impacted point has 
to accelerate from zero, this is, the first term to the left of the equal sign has 
to be zero but for an unidimensional string without bending strength and with 
a mass much small than the mass of the body, its point in contact with the 
material point will have the same velocity when contact each other in 
concordance with the assumption made in this thesis where the contacted 
string point instantaneously achieves a finite velocity after being impacted. 
And here is the contradiction, equation (94) should model a movement where 
the string point has to have a finite velocity but in reality the wire rope cross 
section or axis accelerates from zero when is impacted. Hence the behavior 
of the wire rope under transverse impact from the instant of contact to the 
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time when the impact body moves conjointly with the wire rope section over it 
impacts can not be modeled using an unidimensional string model and leads 
to a contradiction. Thus, here is an explanation: 


Idealizing a wire rope as an string and a solid as a material point 
provokes that after its impact a discontinuity on the tangent of the string. 
Nevertheless the Euler’s general laws of movement can only be used if the 
body considered is continuous and of class C3. And as the idealization made 
here uses an unidimensional string and a material point, it was possible to 
use Newton’s third law because it is only applicable to material points and not 
to continuities. Although Newton never formulated but only enunciated his 
laws, it was Euler who exp,rf)ai^^/»lfef/t^ematically for deformable continuous 


bodies [77], [78], Citin 

/ 

"Every body per6 
ahead, unless it be co 
it." In the generality 
always valid, for a bocy^nay 



UK? OUUJCUL 


rest or of uniform motion straight 
s t state by forces impressed upon 
n mechanics, this axiom is not 
internal or external constraints not 


expressed in terms of a syscam of 'roes. For example, a rigid body subject 
to no applied force spins about some axis through its center of mass; its 
parts, which also are bodies, move in such a way that their centers of mass 
describe circles about that axis. Newton himself did not specify any 
mathematical properties of bodies or forces, and so his intentions must be 
inferred by the reader, and in the course of time different readers have read 
different meanings into his words.” [78] 


In conclusions the foregoing development also means that 
unidimensional idealizations can not be used to model the first instants of the 
impact between bodies (a time that can be called impact period instead to the 
whole interaction) because its not possible to take in count the relative 
velocity between its points in contact and up to them could share a common 
velocity. Furthermore, contact happens between surfaces and not necessarily 
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between the bodies gravity centers. This is why the usual momentum model 
of two spheres impacting can not consider tribological phenomena in and 
around the contact zone. In addition, momentum law used for material points 
only approximates what happens before and after the impact happens and 
therefore is no way to know the strain inside the impacted bodies. 

8.2 . Annex 2. 

Figure 60 shows the results of other 3 numerical models or codes 
developed for many years to study wire ropes under impact (longitudinal and 
transverse) and that were tried by C. Stolle in its thesis to compare the 
effectiveness of his own numerical model and simulation. In the next figure is 
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Time (s) 

Figure 60: Comparison of the mathematical model of the anchorages tension T(KN) applied to the transversely impacted wire 
rope used in the experiment 1 of C. Stolle, jointly with the simulation of C. Stolle and other numerical codes that he used to 
compare his results. [31]. 









